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KAM FOR THE NON-LINEAR BEAM EQUATION 2: 
A NORMAL FORM THEOREM 

L. HAKAN ELIASSON, BENOiT GREBERT, AND SERGEI B. KUKSIN 


Abstract. We prove an abstract KAM theorem adapted to space-multidimensional 
hamiltonian PDEs with regularizing nonlinearities. It applies in particular to 
the singular perturbation problem studied in the first part of this work. 
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1. Introduction 

1.0.1. The phase space. Let A and T be two finite sets in TA* and let Too be an 
infinite subset of lA-. Let T be the disjoint union AU TU Too and consider (C^)^. 
For any subset X of £, consider the projection 

TTX : (C^)^ ^ (C^)^ = {C G (C^)^ : C, = 0 Va ^ X}. 

We can thus write (C^)'^ = (C^)^ x C = (Ca:, C£\x)j and when X is finite 

this gives an injection 

tx : (C^)^ 

whose image is (C^)^. 

Let 7 = ( 71 , 72 ) G and let Ky be the space of sequences ( G (C^)^ such that 

IICII 7 = jElCaPe27da|(«)2,. < 

y ae£ 

- here (a) = max(|a|, 1) and | • | is the standard Hermitian norm on C" associated 
with the standard scalar product (•, •)c"- 

Write Ca = {Pa,qa) and let 

^(C^C') = '^PaQa - QaPa- 
aec 

n is an anti-symmetric bi-linear form which is continuous on 

Yj X F_-y U Y-j X C 

with norm ||11|| = 1. The subspaces are symplectic subspaces of two (com¬ 

plex) dimensions carrying the canonical symplectic structure. 

fl defines as usual (by contraction on the first factor) a bounded bijective operator 

yy9C^f^(C,-)eil%EI 

We shall denote its inverse by 

J : ^ Ky. 

NB. There is another common way to identify YT^ with Yy, the L^-pairing. This 
pairing defines an isomorphism V : YTy —^ Yy such that 

JoV-iC = {( J ~ 0 ^ )Ca:aG/:}. 

The operator is a complex structure compatible with fl which is customarily 

denoted by J, and we shall follow this tradition. This abuse of notation will cause 
no confusion since the two J’s act on different objects: one acts on one-forms and 
the other on vectors, and which is the case will be clear from the context. 

A bounded map A : Fy —>■ Fy, 7 > (0, O),0 is symplectic if, and only if, it extends 
to a bounded map A : F_-y —> F_y, and verifies 

L!(AC,AC')=f^(C,C'), (TYyA'TY.y, 

or, equivalently. A* o J~^ o A = J~^ on Fy and on F_y. If A is bijective, then it is 
symplectic if, and only if. A* o o A = on Fy (see HU)- 

^ Y* denote the Banach space dual of Yj 
^ (7I.72) < (71^72) if, and only if 7J < 71 and 7J < 7^ 
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Let 

= c-^ X (C/27rZ)-^ 

and consider the Banach manifold A'^ x TTc\A^'y whose elements are denoted x = 

(r,6» = [z],w). 0 

We provide this manifold with the metric 

IIx-t'II = inf \\{r,z + 2-Kp,w)-{r',z',w')\\^. 

We provide A-^ x 7r£y_4Ky with the symplectic structure fl. To any C^-function 
f{r,9,w) on (some open set in) A-^ x 'Kc\a^i it associates a vector field Xf = 
— J{df) - the Hamiltonian vector field of / 0 - which in the coordinates {r,0,w) 
takes the form 

V y V J \ da J y j ■ 

1.0.2. An integrable Hamiltonian system in oo many dimensions. In this paper we 
are considering an infinite dimensional Hamiltonian system given by a function 
h{r,w,p) of the form 

(1.1) {r,uj{p)) + ^{w,A{p)w) = {r,uj{p)) + ^{wjr,H{p)wjr) + ^ ^ Xa{pl+ql), 

where Wa = {Pa,qa) and 
r to 

( 1 . 2 ) < Aq ; 2 ? —>■ K., a G Coo 

[ H :V-^ gl{R^ X R-^), *H = H 

are C®*, s, > 1 , functions of p S 2 ?, the unit ball in parametrized by some 
finite subset V of Z'^* . 

The Hamiltonian vector field of h is not on A^ x 'Kc\A^'yi but its Hamiltonian 
system still has a well defined flow with a finite-dimensional invariant torus 

{0} X T-^ X {0,0} 

which is reducible^ i.e. the linearized equation on this torus (is conjugated to a 
system that) does not depend on the angles 0. This linearized equation has infinitely 
many elliptic directions with purely imaginary eigenvalues 

{lAa(p) ■ a tz Aqo} 

and finitely many other directions given by the system 

(jr = JH{p)Cjz. 


^ [z] being the class of z G 

^ there is no agreement as to the sign of the Hamiltonian vectorfield - we’ve used the choice of 
Arnold [T] 
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1.0.3. A perturbation problem. The question here is if this invariant torus for h 
persists under perturbations h + f, and, if so, if the persisted torus is reducible. 

In finite dimension the answer is yes under very general conditions - for the first 
proof in the purely elliptic case see [3], and for a more general case see [ 8 ]. These 
statements say that, under general conditions, the invariant torus persists and 
remains reducible under sufficiently small perturbations for a subset of parameters 
p of large Lebesgue measure. Since the unperturbed problem is linear, parameter 
selection can not be avoided here. 

In infinite dimension the situation is more delicate, and results can only be proven 
under quite severe restrictions on the normal frequencies (i.e. the eigenvalues IAq). 
Such restrictions are fulfilled for many PDE’s in one space dimension - the first 
such result was obtained in m 

For PDE’s in higher space dimension the behavior of the normal frequencies is 
much more complicated and the results are more sparse. A result for the Beam 
equation (which is simpler model than the Schrodinger equation and the Wave 
equation, and frequently considered in works on nonlinear PDE’s) was first obtained 
in [9] and [10]. For other results on PDE’s in higher space dimension see the 
discussion in the first part of this work |5]. 

1.0.4. Conditions on the unperturbed Hamiltonian. The function h we shall consider 
will verify several assumptions. 

Al - spectral asymptotics. There exist constants 0 < c',c < 1 and exponents 
Pi = 2, P 2 > 0, P 3 > 0 such that for all p € V: 


(1.3) 


|A.(p)-|arM<c 


{a)P^ 


d G >Cc 


(1.4) |(A,(p)-A,(p))-(|a|^^ 


\bf^)\< 


< c^cmax( 


(a)/33 ’ {b)A. 


■), 


n, b ^ ^00 , 


\ r ^ai.p') Pl C a (z C.00 

> \ ||(Jif(p))-l||<^; 

N / l(Aa(p) - Ab(p)))| > c' a,b & Coo, \a\^\b\ 

> \\\{\a{p)I-UH{p))-^\\<^ as £00. 

A2 ~ transversality. We refer to section 13.21 for the precise formulation of this 
condition which describes how the eigenvalues vary with the parameter p. 


1.0.5. Conditions on the perturbation. For cr > 0 we let Oj{a,p) be the set 


{x = (r, 0, w) G X TTc\A^i 


r .S 6 » w 
(-, 1—,-)-0 
pap 


< li¬ 


lt is often useful to scale the action variables by p^ and not by p, but in our case 
p will be ~ 1 , and then there is no difference. 

We shall consider perturbations 


f : Oj^{a,p) ^ C, 7 * = (0,m,) > (0, 0), 
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that are real holomorphic up to the boundary (rhb). This means that it gives real 
values to real arguments and extends holomorphically to a neighborhood of the 
closure of (cr, ^). / is clearly also rhb on O^' (a, /r) for any 7' > 7*, and 

Jdf : 0-yi{a, n) — Y-y 

is rhb. But we shall require more: 

R1 - first differential 

Jdf : 0-y'{a,fi) —>■ Ky' 
is rhb for any 7* < 7' < 7. 

This is a natural smoothness condition on the space of holomorphic functions on 
jji), and it implies, in particular, that Jd^f(x) € B(Yy;Yy) for any x € 

That Jd'^f{x) € 13 {Yj;Yy) implies in turn that 

I Jd^/(x)[ea,ef,]| < min(-^, 

{ 0 } {a} 

for any two unit vectors Ca € and Cb € But many Hamiltonian 

PDE’s verify other, and stronger, decay conditions in terms of 

min(|a — 6|, |a + 5 |). 

Such decay conditions do not seem to be naturally related to any smoothness con¬ 
dition of /, but they may be instrumental in the KAM-theory for multidimensional 
PDE’s: see for example [7] where such conditions were used to build a KAM-theory 
for some multidimensional non-linear Schrodinger equations. 

The decay condition needed in this work depends on a parameter 0 < >^ < m* 
and defines a Banach sub-algebra Mfy ^ of B{Y^-, Y^), with norm ^ - its precise 
definition will be given in section [^?T] We shall require: 

R2 - second differential 

Jd^f {a, 

is rhb for any 7* < 7^ < 7s. 

Denote by T-y^i^{cr,p) the space of functions 

/ : C>7. (o'- d-) ‘C- 

real holomorphic up to the boundary, verifying R1 and R2. We provide Ty,i^{o',p) 
with the norm 

sup-y^<7'<7 suPxeOy (ay) \ \Jdf{x)\\Y 

®'^P7.<7'<7 SUP:,eOy (ay) 11 11 7 ' 

making it into a Banach space. Notice that the first two “components” of this norm 
are related to the smoothness of /, while the third “component” imposes a further 
decay condition on Jdf f{x). 
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1.0.6. The normal form theorem. For any a S Coo, let 

[a] ={b^Coo- \b\ = |a|}. 

Theorem. Let h he a Hamiltonian defined by (EH) and verifying Assumptions 
Al-2. 

Let f : (cr,/i) —^ C 6 e real holomorphic and verifying Assumptions Rl-2 with 

7 = ( 7 i,m,) > 7 , = ( 0 ,m*) and 0 <x<m*. 

If e = |/|o-,/i is sufficiently small, then there is a set D' gD with 

Leh{V \V')^0, e^O 

and a C®* mapping 

$ : 0^.(CT/2,/r/2) xV ^ 

real holomorphic and symplectic for each parameter p G D, such that 

(/i +/) o $ =/i'+/'e 

and 

(i) for p gV and ( = r = 0 

drf = def = dcf = dlf = 0; 

(ii) h' : (cr/2,/i/2) x 2? —> C is a 0'^“ -function, real holomorphic for each 

parameter p GV, and 

[d^p{h'{-,p) - < Ce, \j\ < s* - 1 ; 

(iii) h' has the form 

{r,uj'{p)) + ]^{CF,H'{p)Cjr)-G 

+ 9 X] ((P[a]’^[a](p)P[a]) + ('IH, ^[a] (p)<?[a])) 

[a] 

where the matrix is symmetric; 

(iv) for any x G ((t/2 ,/i/2) 

\\dj,i<I>{x,p)-x)\\-y, < Ce, \j\ < s* - 1. 

The constant C depends on ffA, ffT, ffV, d*, s^,m^, >c and h, but not on e. 

We shall give a more precise formulation of this result in Theorem 13.61 and its 
Corollary 13.71 

1.0.7. A singular perturbation problem. We want to apply this theorem to construct 
small-amplitude solutions of the multi-dimensional beam equation on the torus: 

Utt + CAu mu = —g{x, u), u = uit, x), x G T'^*. 

Here g is a, real analytic function satisfying 

g{x, u) = -I- 0{u*). 

Writing it in Fourier components, and introducing action-angle variables for the 
modes in (an arbitrary finite subset) A, the linear part becomes a Hamiltonian 
system with a Hamiltonian h of the form EH, with T void, h satisfies (for all 
m > 0) condition Al, but not condition A2. 
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The way to improve on h is to use a (partial) BirkhofF normal form around u = 0 
in order to extract a piece from the non-linear part which improves on h. This leads 
to a situation where the assumptions A1 and A2 and the size of the perturbation 
are linked - a singular perturbation problem. 

In order to apply the theorem to such a singular situation one needs a careful and 
precise description of how the smallness requirement depends on the (parameters 
determining) assumptions A1 and A2. This is quite a serious complication which is 
carried out in this paper and the precise description of the smallness requirement 
is given in Theorem l3.6l 

This normal form theorem improves on the result in and |10j in two respects. 

• We have imposed no “conservation of momentum” on the perturbation - 
this has the effect that our normal form is not diagonal in the purely elliptic 
directions. In this respect it resembles the normal form obtained for the 
non-linear Schrbdinger equation obtained in [7] and the block diagonal form 
is the same. 

• We have a finite-dimensional, possibly hyperbolic, component, whose treat¬ 
ment requires higher smoothness in the parameters. 

The proof has no real surprises. It is a classical KAM-theorem carried out in a 
complex situation. The main part is, as usual, the solution of the homological equa¬ 
tion with reasonable estimates. The fact that the block structure is not diagonal 
complicates, but this was also studied in for example [7]. The iteration combines a 
finite linear iteration with a “super-quadratic” infinite iteration. This has become 
quite common in KAM and was also used for example in [7]. 

1.0.8. Notation and agreements, (a) = max(|a|, 1). i will denote the complex imag¬ 
inary unit, z is the complex-conjugate of z G C. By (C, we denote the standard 
Hermitian scalar product in C", conjugate-linear in the first variable and linear in 
the second variable. 

All Euclidean spaces, unless otherwise stated, are provided with the Euclidean 
norm denoted by | • |. Eor two subsets X and F of a Euclidean space we denote 

^inf |a;-y| 

rrGA, y^Y 

and 

diam(A) = sup \x — y\. 

x,yGX 

The space of bounded linear operators between two Banach spaces X and Y is 
denoted B{X;Y). Its operator norm will usually be denoted || • || without specifi¬ 
cation of the spaces. Our complex Banach spaces will be complexifications of some 
“natural” real Banach spaces which in general are implicit. An analytic function 
between domains of two complex Banach spaces is called real holomorphic if it gives 
real values to real arguments. 

The sets M, P, £oo,’P, as well as the “starred” constants will be fixed 

in this paper - the dependence of them is usually not indicated. 

Constants depending only on the dimensions on d*,s,,m, and 

on the choice of finite-dimensional norms are regarded as absolute constants. An 
absolute constant only depending on (3 is thus a constant that only depends on 
/3, besides these factor. Arbitrary constants will often be denoted by Ct., ct. and, 
when they occur as an exponent, exp. Their values may change from line to line. 
For example we allow ourselves to write 2Ct. < Ct.. 
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2. Preliminaries. 

2.1. A matrix algebra. The mapping 

( 2 . 1 ) (a, 6 ) I—>■ [a — 6 ] = min(|a — b\, \a + 6 |) 

is a pseudo-metric on lA-, i.e. verifying all the relations of a metric with the only 
exception that [a — b] is = 0 for some a ^ b. This is most easily seen by observing 
that [a-b]= dHausdorff({±a}, {±6}). We have [a - 0] = |a|. 

Define, for any 7 = ( 71 , 72 ) > (0, 0) and > 0, 

e.y,^(a, b) = max([a — 6 ], 1 )'’'^ min((a), (b))’^. 

Lemma 2.1. 

(i) // 7 i ,72 - > 0 , then 

b) < c)e^^^{c, b), Va, b, c, 

if C is sufficiently large (bounded with x)- 
(ii ) If —7 < 7 < 7 , then 

(a, 0 )< (a, 6 )e^^>j(f>, 0), Vu, b 

if C is sufficiently large (bounded with ^ 2 , x). 

Proof, (i). Since [a — 6 ] < [a — c] -I- [c — b] it is sufficient to prove this for 71 = 0. If 
72 = 0 then the statement holds for any C > 1, so it is sufficient to consider 72 > 0 
and, hence 72 = 1. Then we want to prove 

max([a — 6 ], 1 ) min((a), {b)ffi 

< C max([a — c], 1 ) max([c — b], 1) min((c), (b))’^. 
Now max([a — 6 ], 1) < max([a — c], 1) -|- max([c — b], 1), 

max([c — 6 ], 1 ) min((c), {b)ffi > {bffi: 

and 

max([a — c], 1 ) min((c), {b)ffi > min((a), {h)ffi. 

This gives the estimate. 

(ii) Again it suffices to prove this for 71 = 0 and 72 = 1- Then we want to prove 

max(|a|, 1)^^ < C max([a — 6 ], 1) min((a), {b)ffi max(| 6 |, 1)^^. 

The inequality is fulfilled with C > 1 if a or 6 equal 0. Hence we need to prove 

|a|^^ < C'max([a — 6 ], 1 ) min((a), {b)ffi . 

Suppose 72 > 0. If |a| < 2 | 6 | then this holds for any C > 2. If |a| > 2 | 6 | then 
[a —b] > i |a| and the statement holds again for any C >2. 

If instead 72 < 0, then we get the same result with a and b interchanged. 

□ 
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2.1.1. The space We shall consider matrices A ■. C x C ^ gl{2,C), formed 

by 2 X 2-blocs, (each is a 2 x 2-matrix). Define 




= max 


f sup^X^b 


1 SUPfcXla 

Al 


b') 


where the norm on is the matrix operator norm. 

Let denote the space of all matrices A such that \A\^ ^ < oo. Clearly ^ 

is a norm on It follows by well-known results that provided with this 

norm, is a Banach space. 

Transposition - {*A)^ = ^A^ - and C-conjugation - {A)^ = - do not change 

this norm.The identity matrix is in if, and only if, = 0 , and then |/|^,o = C'- 


2.1.2. Matrix multiplication. We define (formally) the matrix product 

{ABt = Y,AlBl 

C 

Notice that complex conjugation, transposition and taking the adjoint behave in 
the usual way under this formal matrix product. 

Proposition 2.2. Let 72 > If A £ o-nd B £ then AB and 

BA £ and 

and \BA\^^^<\A\^,\B\^^^. 

Proof, (i) We have, by Lemma [2dT i') . 

^ |(Ai?)^| e^.^(a, 6 ) < ^ \B^,\ ey,Ua,b) < 

b b,c 

— ^ ^ l^ol \^c I b) 

b,c 

which is < |A|..^ q \B\y This implies in particular the existence of (AB)^. 

The sum over a is shown to be < |A|..^ g |i?|^ ^ in a similar way. The estimate of 
BA is the same. □ 


Hence Aiyp is a Banach algebra, and is a closed ideal in M.yp. 

2.1.3. The space My_^. We define (formally) on Yy (see section IT. 0 . 11 ) 

{A0a=J2A%- 

b 

Proposition 2.3. Let —7 < 7 < 7 . If A G and G then A( G Yj and 

UC\\y<\A\^^JC\\y- 

Proof. Let C = A(. We have 

^ |C7y.o(a,0)2 < E (E l^'l 141 e^,o(a,0))4 

a a b 

Write 

I A^l ICbl e^.o(a, 0) = J X (J |Cb| Cy^oib, 0)) x J 
where _ 

I = Ia,b = flYYIYb)- 
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Since, by Lemma lOT iii . 

j _ e^^o(a, 0 ) ^ ^ 

e^,^(a, 6 )e^,o(&, 0 ) “ 

we get, by Holder, 

^ ICI" eyo(a, 0)2 < /2 J(^ /2 , |C,|2 e^,o(6, O)^) 

a abb 

< E IC^l' eyo(6,0)2 < ^ lai^ eyo(&, 0)2 ^ 4% < 

a,b b a 

This shows that ya exists for all a, and it also proves the estimate. □ 


We have thus, for any —7 < 7 < 7, a continuous embedding of 

"-A AAj^o —y B{Yj] Ty), 

into the space of bounded linear operators on Yx^. Matrix multiplication in 
corresponds to composition of operators. 

For our applications we must consider a larger sub algebra with somewhat weaker 
decay properties. For 7 = (71,72) > 7* = ( 0 ,m*), let 

which we provide with the norm 

11^117,= II^IIb(V.^;K,) + l^l( 7 i. 72 +>f-m.),>< ’ 

This norm makes q into a Banach sub-algebra of B{Yj; Yy) and A 4 y ^ becomes 
a closed ideal in My q. 


2 . 2 . Functions. Let 

and let 


0 < x: < m* 

7 = (71,m*) > 7* = (0,m,). 


2.2.1. The function space 7'y,j<. Consider the space of functions / : O 7 , (cr,/i) —>■ C 
which are real holomorphic up to the boundary (rhb) of Oy^ (cr, /r). This implies 
that for any 7 > 7 * 

/ : Oy{a,fj,) -)■ C 
Jdf : Oy{(T,fJ,) —>■ Y_y 

and 

J(ff : Oy{cr,fi) B{Yy,Y_y) 

are also rhb. We define Ty,:^c{cr, fj.) to be the space of such functions such that for 
any 7 * < 7 ' < 7 , 

Ri — Jdf : Oyi{a,yL) —>■ Yy^ 
i?2 - Jd^f-.Oy'{(T,y)^My,^ 

are rhb. 

We provide 4 ^x((T,/ r) with the norm (11.71) . The higher differentials d^^"^f can 
be estimated by Cauchy estimates on some smaller domain in terms of this norm. 
This norm makes /i) into a Banach space and a Banach algebra with the 

constant function / = 1 as unit. 



KAM FOR THE NON-LINEAR BEAM EQUATION 2 


11 


Remark. The differential forms d^^'^f{x), x G (cr,/i), are canonically identified 
with three bounded linear maps 

> C 


-^0 ®fc-i-2 ~ y~f 

Ri — Y-!' 




fc +2 times 
J 




R2 - ®k ^ 7 ' 
where ® is the symmetric tensor product (over C). 


j 




2.2.2. The function space T-y.jt.v- Let V be the unit ball in R^. We shall consider 
functions 

/ : O-y. {(7,fi)xV^C 

which are of class C®*. We say that / G h) if, and only if, 

f 

for any p G 7) and any |j| < s*. We provide this space by the norm 

d^f 

|/|<T,ii = max sup I-—— ("i 

This norm makes /i) into a Banach space and a Banach algebra. 

2.2.3. Jets of functions. For any function / G 7fy,>f,r>((7,/i) we define its jet /^(x), 
X = (r, 6, w), as the following Taylor polynomial of / at r = 0 and w = 0 

(2.2) /(O, 0, 0) + drfiO, 0, 0)[r] + d„/(0, 0, 0) [w] + ^4/(0, 0, 0) [w, w] 

Functions of the form will be called jet-functions. 

Proposition 2.4. Let f G 71),>r,-D(CT, ^). Then G 71(.,:«r,r>(cr,/r) and 

l/^l <C\f\<y,ii . 

7 ,>tr,X> 7 ,x,X> 

C is an absolute constant. 

Proof. The first part follows by general arguments. Look for example on 
g{x) = dlf op{x)[w,w], x={r,0,w), 

where p{x) is the projection onto (0,0,0). This function g is real holomorphic up 
to the boundary (rhb) on (cr,/i), being a composition of such functions. Its 
sup-norm is obtained by a Cauchy estimate of /: 

\\dlfiPix))\\^,Y V -C) ll^^lly ^ Ct.^ sup |/(y)| ||w||^ < Ct. sup \f{y)\. 
Since Jdg{x)[-] equals 

{Jddlf op{x)[w, ti;]) [dp[-]] + 2{Jdl,f o p{x)[w]) [•], 

and 

Jdlf 

and 

Jddlf = Jd^df : Oy{a,p) B{Y^>,Yy]Y^>) 
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are rhb, it follows that dg verifies R1 and is rhb. The norm || J(i( 7 (x)||^, is less than 
pdldf{p{x))\\i^^y ^ Y ||w||^, + 2 \\jdlf{p{x))P^y llwll^,, 


which is < Ct. \\Jd,f{y)\\^, - this follows by Cauchy estimates of deriva¬ 

tives of Jdf. 

Since Jd‘^g{x)[,-] equals 

{Jd'^dlf op{x)[w, w]) [dp[-], dp[-]] + 2J{ddl,f op{x)[w]) [■,dp[-] + 2Jdlf op{x)[-,-], 


and 

Jddlf = Jdl,df : -A ^(T^/;7Wy 

and 

Jd^dlf = Jdldp : Oy(a,Ai) ^ 

are rhb, it follows that Jdg^ verifies R2 and is rhb. The norm || ^ is less 

than 


Jd^d fipix))\ 


B(Yy,Yy,M’’ 


kli: 


\jd^,d'^f{p{x))P^y j llwll. 


-f2||jdV(a;)||^,^^, 

which is < Ct.supygg, |||| / ” this follows by a Cauchy estimate of 

Jd^f. 

The derivatives with respect to p are treated alike. □ 


2.3. Flows. 


2.3.1. Poisson brackets. The Poisson bracket {f,g} of two C^-functions / and g is 
(formally) defined by 

n{Jdf, Jdg) = df[Jdg\ = -dg[Jdf] 

If one of the two functions verify condition Rl, this product is well-defined. More¬ 
over, if both /a nd g are jet-function, then {f,g} is also a jet-function. 


Proposition 2.5. Let f^g £ m), and let a' < a and p' < p <1. Then 

(i) and 




for 




1 


(cr-cr') {p-p')^' 
(ii) the n-fold Poisson bracket Pff f £ 7i),>c,xi(cr, p) and 

where Pgf = {g,f}. 

C is an absolute constant. 
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Proof, (i) We must first consider the function h = Q{Jdg, Jdf) on O 7 . (cr, /r) Since 
Jdg, Jdf : 0.y^{a,g,) —?► liy, are real holomorphic up to the boundary (rhb), it 
follows that h : O-y, (cr, /r) —>■ C is rhb, and 

\hix)\ < ||Jd5(a:)ll^. \\Jdfix)\\^^ . 

The vector Jdh{x) is a sum of 

Jil{Jd^g{x),Jdf{x)) = Jd‘^g{x)[Jdf{x)] 

and another term with g and / interchanged. Since Jd^g : 0.yi{a, g,) —>■ B{Y.yi] Yy) 
and Jdg, Jdf : Oji{a,g) —>■ Ty' are rhb, it follows that Jdh verifies R1 and is rhb. 
Moreover 

\\jd^g{x)[Jdf{x),-]\\^, < II \\Jdf{x)\\y 

and, by definition of ^ , 

\\Jd^g{x)\\B^Y_^,^Yy) ^ Pd^9ix)\\Y,o- 

The operator Jdfh(x) = d{Jdh){x) is a sum of 

Jd^g{x)[Jdf{x)\ 

and 

Jd^g{x)[Jd^ f{x)\ 

and two other terms with g and / interchanged. 

Since Jdf'g : —>■ B{Y^r,Mf^, and Jdf : 0^'{a,g) —>■ Ty' are rhb, it 

follows that the first function O^' (ct, g) —S(Ky'; Aly is rhb. It can be estimated 
on a smaller domain using a Cauchy estimate for JdJgix) 

The second term is treated differently. Since 

Jd^f, Jd^g : 0^'{a, g) 

are rhb, and since, by Proposition [221 taking products is a bounded bi-linear maps 
with norm < 1, it follows that the second function 0.y'{a,g) —>■ ^ is rhb and 

\\jd^g{x)[Jd^f{x)]\\_^,^^ < ||Jd 2 g(a;)||^,^^|| JdV(a:)||y^^- 

The derivatives with respect to p are treated alike. 

(ii) That P^f G g) follows from (ii), but the estimate does not follow 

from the estimate in (ii). The estimate follows instead from Cauchy estimates of 
n-fold product Pgf- □ 


Remark 2.6. The proof shows that the assumptions can be relaxed when g is a 
jet function: it suffices then to assume that g G 71 ),0 ,vi<7,g) and g - g(-, 0 , - jl £ 
g). Then {g,f} will still be in 71 ),j<,t)((T,/ r) but with the bound 


7,x,I5 


+ l5-ff(u0, 

7,>f,27 


)l/l 


< 7 ,fl 


To see this it is enough to consider a jet-function g which does not depend on 
9. The only difference with respect to case (i) is for the second differential. The 
second term is fine since, by Proposition 12.21 is a two-sided ideal in Aiy q 

and 


\jd'^g{x)[Jd'^f{x)]\\^,^^ < \\Jd'^g{x 


C'.o 


|7d"/(a 


17 ,>c 


^ g{-, 0, •) this is the 0:th Fourier coefRcent of the function 9 g{-, 9, •) 
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For the first term we must consider Jcfig{x)[Jdf (x)] which, a priori, takes its values 
in MyQ and not in But since g is a jet-function independent of 9 this term 

is = 0 . 


2.3.2. Hamiltonian flows. The Hamiltonian vector field of a C^-function g on (some 
open set in) Ky is —Jdg. Without further assumptions it is an element in FT-y, but 
if g S 71),then it is an element in Yy and has a well-defined local flow 


Proposition 2.7. Let g G ^), and let a' < a and /i' < /i < 1. If 


\g\ < -r; min(cr - a', g - g'), 

7,>f,X7 O 


then 


(i) the Hamiltonian flow map is, for all |f| < 1 and all 7» < 7 ^ < 7 , 

a C®* -map 

Oy {o', g') X V ^ Oy! (tr, g) 

which is real holomorphic and symplectic for any fixed p GD . 

Moreover, 


\dl{^\x,p) -a;)|| < C\g\^,t^ 


-y.x.T) 


and 


\\d^,{d^\x)-I)\\.<C\g\.,n 


.V 


for any x G Oyfla', g'), 7 * < 7 ' < 7 , and 0 < |jj < s*. 
(ii) / o e 7).,x(cr', g', V) for \t\ < 1 and 


\fo<i>llpy <c\f\.,n. 


S|o- .Ai 
'y.x.'D 




C is an absolute constant. 


Proof. It follows by general arguments that ^ : C/ —>■ Oy{a, g) is real holomor¬ 

phic in (t, C) G [/ C C X Oy{a, g) and depends smoothly on any smooth parameter 
in the vector field. Clearly, for \t\ < 1 and x G Oy(cr', g') 

||$‘(a:,p)-a;|| < sup \\Jdg{x)\\ <\g\<y.,, 

^ xe0.y{cr,tJ.) 7,0.15 

as long as <l5*(a^) stays in the domain Oy{a,g). It follows by classical arguments 
that this is the case if 

\g\cr,ij. < ct.rmifla — a', g — g'). 

7,0,T> 

The differential We have 

= -Jd^g{^\x))d^\x) = B{t)d^\x), 

where B{t) G My re-writing this equation in the integral form d^*{x) = 

Id-I- fg B(s)d^^(x)ds and iterating this relation, we get that dd>*(x) — Id = B°°{t) 
with 
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We get, by Proposition 12.21 that — Id G and, for |t| < 1 , 



k>l 

In particular, A = d^*(x) is a bounded bijective operator on Vy. Since Jd^g is 
a Hamiltonian vector field we clearly have that 


Q(A(,Aa = n(c,C), vc,CeYy, 

so A is symplective. 

Parameter dependence. For |j| = 1, we have 



B{t)Z{t)+A{t). 


Since 



it follows by classical arguments, using Gronwall, that 


\\Zm^,<Ct.\g\.,^. \t\. 


The higher order derivatives (with respect tp p) of $*(a;,p), and the derivatives of 
d^*{x,p) are treated in the same way. 

The same argument applies to any 7 » < 7 ' < 7 . 

Since 



n>0 

(ii) is a consequence of Proposition 12.bl iil. 


□ 


3. Normal Form Hamiltonians and the KAM theorem 


3.1. Block decomposition, normal form matrices. In this subsection we recall 
two notions introduced in [7] for the nonlinear Schrodinger equation. They are 
essential to overcome the problems of small divisors in a multidimensional context. 
Since the structure of the spectrum for the beam equation, {\/|a|^ + rn, a G 
is similar to that for the NTS equation, {|ap + a G then to study the 

beam equation we will use tools, similar to those used to study the NTS equation. 

Block decomposition: For any A G N U {c»} we define an equivalence relation 
on Z'^*, generated by the pre-equivalence relation 



(see (j2.1|) '). Let [oJa denote the equivalence class of a - the block of a. For further 
references we note that 


(3.1) 


|a| = l&l and [oJa 7 ^ [&]a => [a — 6] > A 


The crucial fact is that the blocks have a finite maximal “diameter” 


d/\ = max [a — 6 ] 


which do not depend on a but only on A. 
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Proposition 3.1. 

(3.2) dA<CA^^. 

The constant C only depends on d*. 

Proof. In [7] it was considered the equivalence relation on generated by the 
pre-equivalence 

atv b if |a| = |6| and |a — 6| < A. 

Denote by [a]^ and the corresponding equivalence class and its diameter (with 
respect to the usual distance). Since a ^ 6 if and only if a ~ 6 or a ~ —b, then 

(3.3) [a]A = NaU-Ha, 

provided that the union in the r.h.s. is disjoint. It is proved in [7] that d^ < Da ='■ 
CA' ' . Accordingly, if |a| > Da, then the union above is disjoint, (IQ) holds 
and diameter of [a] a satisfies (13.2p . If |a| < Da, then [a] a is contained in a sphere 
of radius < Da- So the block’s diameter is at most 2D a- This proves ()3.2I1 if we 
replace there Cd, by 2Cd., ■ □ 


If A = oo then the block of a is the sphere {b : |6| = |a|}. Each block decom¬ 
position is a sub-decomposition of the trivial decomposition formed by the spheres 
{|a| = const}. 


Normal form matrices. Let £a be the decomposition of C = J- U Coo into the 
subsets 

[oJa n Coo a € Coo c Z'’** 

T aG D. 


Ha = 


Remark 3.2. Now the diameter of each block [oJa is bounded by 

if moreover we let C > ffD. 


(d.+l)! 

dA < CA^— 


On the space of 2 x 2 complex matrices we introduce a projection 
n : Mat(2 x 2,C) ^C/-bCJ, 

orthogonal with respect to the Hilbert-Schnhdt scalar product. Note that C/ -I- C J 
is the space of matrices, commuting with the symplectic matrix J. 

Definition 3.3. We say that a matrix A : C x C ^ Mat(2 x 2,C) is on normal 
form with respect to A, AgNU{oo}, and write A G MiF a, if 

(i) A is real valued, 

(ii) A is symmetric, i.e. A^ = 

(hi) A is block diagonal over £a, *.e. A^ = 0 if [oJa ^ Ha; 

(iv) A satisfies IIA^ = A'^ for all a,b G Coo- 

Any quadratic form q(ir) = ^{w, Aw), w = {p, q), can be written as 
i(p, Aiip) -h (p, Ai 2(?) -b i((7, A22q) + ^{wjr,H[p)wjr) 
where An, A 22 and H are real symmetric matrices and A 12 is a real matrix. 
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We now pass from the real variables Wa = (jpa,qa) to the complex variables 
Za = (^a, ?7a) hy w = Uz defined through 


(3.4) 




for a € Coo^ and acting like the identity on Then we have 


^{Uz) = + itQv) + 

where 

P = -^ ((^11 ~ ^22) — i(^i 2 + *^12)) 

and 

< 5=2 ((^11 + ^22) + i(^i 2 — *^12)) ■ 

Hence P is a complex symmetric matrix and Q is a Hermitian matrix. If A is on 
normal form, then P = 0. 

Notice that this change of variables is not symplectic but 




iJoo 0 \ 

0 J:f ]' 


3.2. The unperturbed Hamiltonian. Let h he a function as in (HH), i.e. 

(3.5) h{r,w,p) = (r,w(p)) + ]^{w,A{p)w), 

where 

(w,A{p)w) = {wj7,H{p)wjr) + ^{{poo,Q{p)p oo) T {QoO 5 Q{p)qoo)) 

and 

Q{p) = diag{Aa(p) : a € Coo}- 
Assume w, Xa, H verify (O. We shall denote by 

X = |Vpa;|c>!.-i(x>) + sup |VpAo|c«.”i(T>) + pH\\c‘*-'^CD) 

aec 

since this quantity will play an important role in our analysis. 


3.2.1. Assumption A1 - spectral asymptotic. 

There exist constants 0 < c, c' < 1 and exponents /3i = 2, /32 > 0, /Sa > 0 such 
that for all p G P the relations (HSl), (HI, dfll) and (11.31) hold. 


3.2.2. Assumption A2 - transversality. Let [a] = [a]oo so that [a] equals {b : | 6 | = 
|a|} when a € Coo and equals P when a G T . Denote by Q[a] the restriction of the 
matrix Q to [a] x [a] and let ( 5 [ 0 ] = 0. Let also JP(p)[ 0 ] = 0 and m = 2#P. 

There exists a 1 > i5o > 0 such that for all C®*-functions 
oj' : D ^ R", \(xi' — w|cs.(x)) < Aq, 

the following hold for each fc G Z” \ 0 : 
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(i) for any a,b G Coo U { 0 } let 

L{p) : X ^ {k,u;'ip))X + Q[a](p)^ ± XQ^.y, 
then either L{p) is 5o-invertible for all p G V , i.e. 

oo 

or there exists a unit vector 3 such that 

1(11,332.(^)11)1 > ( 5 o, ypGiE 

and for any unit-vector v in the domain of L{p)] 

{a) let 

L{p, A) : X hA (fc, Lo'{p))X AX iXJH{p) 

and 

P{p, A) = det L{p, A) : 

then either L{p, \a{p)) is (5o-invertible for all p gV and a € [a] 00 , or there 
exists a unit vector 3 such that 

\diP{p,K{p)) + dxP{p,Xa{p)){v,d^Q{p)v)\ > Ao WH', 

for all p € 2i>, a G [a]oo and for any unit-vector v G (C^)[“l; 

(ill) for any a,b G PU {0} let 

L(p) : X ^ (k,u;'(p))X - iJII(p)ia]X + iXJH(p)ii,y 

then either L(p) is Ao-iRvertible for all p G P, or there exists a unit vector 
3 and an integer 1 < j < s* such that 

\dl deti(p)| > (5o ||i.(p)||c.7D), 'ipGV. 


Remark 3.4. The dichotomy in A2 is imposed not only on uj but also on C^- 
perturbations of uj' , because, in general, the dichotomy for uj' does not imply that 
for perturbations. 

If, however, any C* perturbation uj' of uj can be written as uj' = uj o f for some 
diffeomorphism / = id + 0{5{)) - this is for example the case when uj{p) = p, - then 
the dichotomy on uj implies a dichotomy on C'**-perturbations. 


3.3. Normal form Hamiltonians. Consider now the function (13.51) defined on 
the set T). This function will be fixed throughout this paper and we shall denote it 
and its “components” by 


lup, ^up: Aup, Qup^ Pu 


Remark 3.5. The essential properties of /lup are given by the constants 

X,c',c,/3 = (/3 i,/32,/33),^o- 

These will be fixed now, once and for all. All estimates will depend on h^p only 
through these constants. Since it will be important for our analysis of the Beam 
equation we shall track this dependence with respect to c', Sq, x- In order to simplify 
the estimates a little we shall assume that 


(3.6) 


0 < c' < Jo < X ^ c. 


® O3 denotes here the directional derivative in the direction 3 € 
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We shall consider functions of the form 

(3.7) h{r, w, p) = {u}{p), r) + i(w;, A{p)w) 
which satisfies 

Hypothesis w: oj is of class C®* on V and 

(3.8) \uj — Wuplcfl,(X)) < (5. 


Hypothesis B: A — Aup : T> —?> A4q ^ is of class C®*, A(p) is on normal form 
G for all p G 27 and 

(3.9) \\dj,{A{p) - AM)la]\\< 

for |j| < s*, a G £ and p G 27. Q We also require that 

(3.10) x>0. 

A function verifying these assumptions is said to be on normal form, and we 
shall denote this by 

Since the unperturbed Hamiltonian hup will be fixed in this paper we shall often 
suppress is in this notation writing simply h G MT5). 

3.4. The normal form theorem. In this section we state an abstract KAM result 
for perturbations of normal form Hamiltonians by a function in h), 0 < 

(T, /i < 1. Let 

7 = ( 71 , m*) > 7 * = ( 0 , m*) and > 0 . 


Theorem 3.6. There exist positive constants C, a, exp and expg such that, for 
any h G (A, (5) and for any f G Tf,^,vicr, p), 

£ = 1/^1 and f = |/| o-,p 




7,>r,'D 


*/ 

and 

(3.11) 


a < —r' 
- 2C 


e(logi)-P < l( max(7, \dA) ^-exp.^.exp3^,^g 
£ ap ' X + s 


then there exist a set 27' = 'D'{h,f) C 27, 

Leb(I) \ BO < C( log ' p, 

e: (7/i (^0 ^0 

and a C^* mapping 

$ : O^S^I2,pI2) X 27 -)► 0^,{a,p), 
real holomorphic and symplectic for each parameter p G 27, such that 

(h + /) o $ = h' + /' 

and 


^ here it is important that 11 11 is the matrix operator norm 
® remember the convention (13.611 
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(i) for p gV' and ( = r = 0 

drf = dsf = dcf = dlf = 0; 

(ii) h' G MJ-S' < yj o.iT'd 

\d' — ^|cr/2,/i/2 < C'; 

7* 

(iii) for any x G 0-y^(u 12, pl2) and |j| < s* 

(iv) if p = (0, p 2 , ■ ■ ■, Pp) and f"^{-,p) = 0 for all p, then h' = h and ■) = x 
for all p. 

C is an absolut constant that only depends on ft, k, c and supp |iu;|. The exponent 
exp' is an absolute constant that only depends on f and >c. The exponent expg only 
depends on s*. The exponent a' is a positive constant only depending on 

The condition on $ and h' — h may look bad but it is not. 


Corollary 3.7. Under the assumption of the theorem, let e* be the largest positive 
number such that (13.111) holds. Then, for any p GT) and |j| < s, — 1, 


{izY 


{my 


9pih'{-,p) - h{-,p))\,^/2,p/2 

7* 



(7,fi ; 

7,>f,X> 


\\di{^{x,p) - x) 


|a^(d$(x,r) - J)||, 


<^\P 

e* 


7,x,D 


for any x G 0.y^{al2, pl2). 

The constant C is an absolute constant that also depends on f). 


Proof. Let us denote p here by p. If |/^| o-./i < then we can apply the theorem 

to ef for any |e| < 1. Let now p = {e,p) and consider h and hup as functions 
depending on this new parameter p - they will still verify the assumptions of the 
theorem, which will provide us with a mapping il> with a C®* dependence in p = (e, p) 
and equal to the identity when e = 0. The bound on the derivative now implies 
that 

(j 

||<i)(a;, £, p) — a;||g,, <Ce< —1/^| o-./i 

for any x G Oj, {<xl2., p/2). The same estimate holds for all derivatives with respect 
to p up to order s* — 1 . 

The argument for h' — h is the same. □ 


We can take 5 = 0 here, in which case h equals the unperturbed Hamiltonian 
hup - this is the case described in theorem of the Introduction. 

Remark 3.8. The assumption (13.111) on e involves many constants and parameters. 
Then (13.111) takes the form 

e(log-rP < 

e - ^x + 7 
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where C depends on ( 77 , ct, /x, A. If we assume that 

X = O((5o"^) and c' = (7((5(5+^) 
for some H > 0, then assumption (13.111) reduces to 

e{log-y^P < c"5i+^+2Ne’'P3^ 

which implies 

^ ^ expg 

when So is sufficiently small. 

Actually in paper Theorem 13.61 this is used in this context. 


4. Small Divisors 

For a mapping L : 2? —>■ (dim, R) define, for any k > 0, 


Let 


S(L,K) = {peV: \\L-^\\ > 


h{r,w,p) = {r,uj{p)) + ^{w,A{p)w) 


be a normal form Hamiltonian in AfT:^c{A, S). Recall the convention in Remark 13. 5 
and assume > 0 and 




(4.1) 

where C is to be determined. 

Lemma 4.1. Let 

Lk = {k,uj{p)). 

There exists a constant C such that if holds, then 

,X + S K 


Leb( U < (77V'="P- 


0<|fe|<IV 


oo Co 


and 

dist(P \ E(Lfe, xc), E(Lfe, ^)) > J 

for any k > 0. 

The exponent exp only depends on ffA. C is an absolut constant. 

Proof. Since <5 < So, using Assumption A2(i), with a = b = 9, we have, for each 
k y 0, either that 

|(a;(p), k)\ > So > K 'ip GV 

or that 

9}{^^{p),k) >So ipeV 

(for some suitable choice of a unit vector 3 ). The first case implies E(Lfc,K) = 0. 
The second case implies that E(Lfe, k) has Lebesgue measure 

So So 


® this is assumed to be fulfilled if ^Lj (§) = 0 
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Summing up over all 0 < |fc| < TV gives the first statement. The second statement 
follows from the mean value theorem and the bound 

\VpLk{p)\<N{x + S). 

□ 


Lemma 4.2. Let 

^fe,[a] = ~ [a]' 

There exists a eonstant C such that if dHH) holds, then, 

Leb( U <CiV-P(^)(^)i 

0<|fe|<A ^0 

[a] 

and 

dist ( 2 ? \ , k), S(Z/fc [q] I ■^)) > 'q ]\[^-^ _|_ ’ 

for any k > 0. 

exp only depends on ^ and ffA. C is an absolut constant that depends on c 
and sup 25 |a;|. 

Proof. Consider first a £ Coc- Then decouples into to a sum of two Hermitian 
operators of the form 

{k, Uj)I + Q[a] 

- denoted L = - where Q^a] is the restriction of Q to [a] x [a]. Since L is 

Hermitian: 

\\L{p)-^\\ < max||((fc,a;(p)) + A(p))"^||, 

where the maximum is taken over all eigenvalues A(p) of Q{p)', 

• for any po G 

9i{v{p)L{p),v{p))/p=p„ = dj{v{po)L{p),v{po))/p=po 

for any eigenvector v{p) of L{p) (associated to an eigenvalue which is in 
the direction 3 ). 

If we let 

Lup — {k, io)I T (Qup)[o]; 

where Qup comes from the unperturbed Hamiltinonian, then it follows, from (13.911 
and mi, that 

\\L — -^'upllcs, (x>) — ^ — ~2 

and, hence, 

dllausdorff (o(L), Cr(Z/up)) < (5. 

If now Lup is do-invertible, then this implies that L is ^-invertible. 

Otherwise 

di{{k,uj{p)) + X{p))/p=p„ = {v{p),d^L{p)v{p)) / p^p„ = 

= {v{p),diLup{p)v{p))/p^p„ +0{6), 
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where v(p) is a unit eigenvector of L{p) associated with the eigenvalue X{p). If L^p 
is not (5o-invertible, then, by assumption A2(i) there exists a unit vector 3 such that 

\di{{k,uj{p)) + X{p)) \ / >So-5>^. 

Hence, the Lebesgue measure of is < “ recall that, by Remark 

13.21 the operator is of dimension < . (This argument is valid if X{p) is in the 

direction 3 which can always be assumed when Q is analytic in p. The non-analytic 
case follows by analytical approximation.) 

Since \(k,uj{p))\ < \k\ < N, it follows, by (I1.3L that 

\{k,uj{p)) + A(p)| > |Aa(p)| - (5 - Ct. |fc| > \af^ - c{a)~^^ - 6 - Ct. \k\ 

for some appropriate a £ [a]. Hence, S(Lfc_[a], k) = 0 for |a| > iV^. 

1 

Summing up over all 0 < |fc| < and all |a| < gives the hrst estimate. 

Consider now a £ Then L{p) = ((fc,u;)/ — iJid) and it follows, by (13.91) and 
(j4.1[) . that 

11^ - -^upllc.. <S< i(5o, 

where L^p{p) = (^{k,ui)I — iJH^p^. If now L^p is do-invertible, then L will be 
^-invertible. 

Otherwise, 

|detL - detLupl^j < Ct.d( ULupH^,- -I- d)"* ^ 

and, by assumption A2(iii), there exists a unit vector 3 and an integer 1 < j < s* 
such that 

\dl detLup(p)| > do \\L^p\\^~f ^-^, Vp G T>. 

This implies that |Tup|cj ^ ct.do and, hence, 

IdetL-detLupIcj < Ct.d ||Lup||(()“^. 


Thus 

and do — Ct.d > 4r. 


dl detT(p)| > (do - Ct.d) , Vp G P, 


— 2 


Then, by Lemma FA. 11 


detL(p) 


1 > K, 

\m — l — ’ 


\\L\\c^ 

outside a set E' of Lebesgue measure 


< Ct.i 


do 




Hence, by Cramer’s rule, 

LebE(L.e) < (A)I < c,A2L±i)(A)!, 

do do do do 

Summing up over all |fe| < N gives the first estimate. 


The second estimate follows from the mean value theorem and the bound 


VpLfc,[a](p)| < N{x + S). 

□ 
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Lemma 4.3. Let 

Lk,[a],[b] = ((^,w)/- iadjA)[a]- 

There exists a eonstant C such that if dHH) holds, then, 

0<W<N 
\a — b\<A' 

and 

dist(P \ S(Z/;j [£,] [f,], k ), S(Z/^ [£,] [f,]! ■^)) ^ N(^x “ 1 “ ’ 

for any n > 0. 

The exponent exp only depends on ^ and ffA. The exponent a is a positive 
constant only depending on Si,, C is an absolut constant that depends on c 

and sup 25 |a;|. 


Proof. Consider first a,b € F. This case is treated as the operator L[p) = ((fc, uj)I— 
iJH^ in the previous lemma. 

Consider then a € Coo and b € F, so that 


Lk,[a]{p)-X^{k,uj{p))X + Q[,]{p)X + XiJH{p). 


Let 

L{p, A) :X^ {k, uj{p))X + AX + iXJH{p) 

and 

P{p, A) = det L{p, A). 

Since [„](«:) is “partially” Hermitian, 

l|ifc.[o](p)"^|| < max||L(p, A(p))“^|| , 


where the maximum is taken over all eigenvalues A(p) of Q{p). 

If we let 

Cupip, A) : X I—>■ {k, ui{p))X + XX + XiJHup{p), 
then it follows, from (13.91) and (|4.1I) . that 

IIL — Lupllc®, (x>) — ^ — ~2 ’ 

If Lup is ^g-invertible, then this implies that L is ^-invertible. 

Otherwise 

^P{p,X{p)) = d^P{p,X{p)) + dxP{p,X{p)){v{p),d.^Q{p)v{p)) = 

= a3Pup(p,Aa(p))+5APup(p,Aa(p))(i;(p),a3Qup(p>(p))+0(^||L||^x7^)). 
By Assumption A2(ii) there exists a unit vector 3 such that 




- 2 


m —1 


Hence, the Lebesgue measure of is ^ d‘^~ recall that, by Remark 

13.21 the operator is of dimension < . (This argument is valid if X{p) is in the 

direction 3 which can always be assumed when Q is analytic in p. The non-analytic 
case follows by analytical approximation.) 
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Let X{p) be an eigenvalue of Q[a]- Since \{k,u}{p))\ < \k\ < N, it follows, by 
(ll.3[l . that 

\{k,u}{p)) + A(p)| > |Aa(/o)| - ^ - Ct. |fc| > |a|^^ - c{a)~^'^ - 6 - Ct.\k\ 

for some appropriate a € [a]. Hence, , k) = 0 for |a| > (N)'^. 

1 

Summing up over all 0 < |fc| < and all |a| < (7V)^ gives the first estimate. 


Consider finally a,b € Loo- Then Lfcja],[&] decouples into a sum of four Hermitian 
operators of the forms 

Lk,[a\,[b]{p) '■ X I—>■ {k,Uj)X -I- (5[a]X -I- 

and 

Lk,la],lb]{p) ■ X l-A {k,Uj)X + Q[a\X - XQli,]- 
The first one is treated exactly as the operator X i-A {k, uj)X+Q[a\X in the previous 
lemma, so let us concentrate on the second one. It follows as in the previous lemma 
that the Lebesgue measure of k') is < ^'a* ^ ~ recall that, by Remark 

13.21 the operator is of dimension < . 

The problem now is the measure estimate of |J k) since, a priori, there 

may be infinitely many that are non-void. We can assume without 

restriction that |a| < \b\. Since \{k, uj{p))\ < Ct. \k\ < Ct.fV, it is enough to consider 
|6| — |a| < Ct.A^ (because /3i > 1). 

Since /3i = 2, |a|^^ — \bf^ is an integer, and outside a set I](2 k') of Lebesgue 
measure 

< Ct.A^-^ 
oo 

we have 

\{k,uj{p)) -h \af^ - \bf^ I > 2k'. 

Then, by (HH), 

\{k,uj{p)) +a{p)-j3{p)\> 2 k!-2-^-2^^ 

\a) 

for any Q;(p) and /3(p), eigenvalues of Q[a]{p) and Q[b]{p), respectively. Now this is 
> k' unless 

|a|<Ct.min((4)i,(4)^)=:^- 

\ Ki Ki / 

Hence, if k' > k, then 

[J S(Lfc_[a].[b], k) C S(2k') U [J E(LfeJa]J;,],K). 

[a],[b],fe |a|,|6|<M-|-Ct.Af,fc 

This set has measure 

^ i5o V k' / (5o 

By an appropriate choice of k' S [k,(5o], this becomes 

< Ct.fV®"P(-^)“ 

Oo 

for some a > 0 depending on 

□ 
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5. Homological equation 
Let /i be a normal form Hamiltonian (j3.7L 

h{r, w, p) = {u}{p), r) + A{p)w) £ NT5) 
and assume >^ > 0 and 

(5.1) 5 < ^c', 

where C is to be determined. Let 


7 = (7,m*) > 7* = (0,m,). 

Remark 5.1. Notice the abuse of notations here. It will be clear from the context 
when 7 is a two-vector, like in H-H^ and when it is a scalar, like in e^'^. 

Let / £ /i). In this section we shall construct a jet-function S that 

solves the non-linear homological equation 

(5.2) {h^s} + {f-f,S}^ + f^ = 0 

as good as possible - the reason for this will be explained in the beginning of the 
next section. In order to do this we shall start by analyzing the homological equation 

(5.3) {h,S} + f = {). 

We shall solve this equation modulo some “cokernel” and modulo an “error”. 


5.1. Three components of the homological equation. Let us write 

f{e,r,w) = frir,0) -b {U{0),w) -b ^{Uw{9)w,w) 

and recall that, by Proposition 12.41 £ T-f,}<,v{cr, p)- Let 

S{9, r, w) = Srir, 9) -b {Su,{9), w) + ]^{Sy,u,{9)w, w), 

where fr and Sr are affine functions in r - here we have not indicated the dependence 
on p. 

Then the Poisson bracket {h, equals 

- {d^Sr{r,9) + {du,Sw{9),w) + ^{dcjSruw{9),w)-\- 

+ {AJSru{9),w) -b -{AJSww{9)w,w) - -{Syjw{9)JAw,w) 

where 3^ denotes the derivative of the angles 9 in direction w. Accordingly the 
homological equation (15.31) decomposes into three linear equations: 

r d^Sr{r,e) = fr{r,9), 

{ d^S^{9)-AJS^{9) = U{9), 
dojSyjW {9)-AJS^n,{e) + Sww {e)jA 

- fww (9). 
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5.2. The first equation. 

Lemma 5.2. There exists constant C such that if (EH) holds, then, for any N > 1 
and K > 0, there exists a closed set = 'Di{h, k, N) C T, satisfying 

LehiV \Vi)< 

Oo Oo 

and there exist C®* functions Sr and Rr on x x 2? —>■ C, real holomorphic 
in r, 9, such that for all p € "Di 

(5.4) d^(p)Sr{r,9,p) = fr{r,e,p) - fr{r,0,pj3- RriO,p) 

and for all {r,9,p) € x T^, x T), |r| < p, a' < a, and \j\ < s* 


(5.5) 


\diSr{r,e,p)\ <C 


1 /ArX +j-Ti 

„-{<7-<7')N 


(5.6) \dj,Rrir,0,p)\<C- - —\f^\a,^ 

\(j — a ) 

The exponent exp only depends on n = ft^A, and C is an absolut constant. 
Proof. Written in Fourier components the equation (15.4|) then becomes, for k € 
Lk{p)S{k) =: {k,uj{p))S{k) = -KF{k) - R{k)) 

where we have written S, F and R for Sr, (fr — fr) and Rr respectively. Therefore 
(15.4[) has the (formal) solution 


with 


and 


5'(r,6»,/9) = and R{r,e, p) F{r,k, p)e^^'^’'^'> 

o(r k o) = i -^k{p)~'^iF{r,k,p) if 0 < |fc| < iV 
S{r,k,p) I Q 

Sip k o') = { 1*1 ^ ^ 

R[r,lt,p) I Q 


By Lemma l4T 


\\iLkip))-^\\< 


for all p outside some set E(Lfc, k) such that 


^{V \ S(Lfc, k), S(Lfe, ^)) > 


and 


Vi^V\ U E(Lfe,K) 

0<|fe|<Ar 

fulfills the estimate of the lemma. 

For p ^ E(Lfc, D we get 

\S{r,k,p)\ < Ct.-\F{r,k,p)\ . 

K 

Differentiating the formula for S{r,k,p) once we obtain 

dj,Sir,k,p) = (^ - -r^dj,F{r,k,p) + {dj,uj,k)F{r,k,p)^ 


10 


/r(r, 0,p) is the 0:th Fourier coefhcent, or the mean value, of the function 0 i—>• fr{r,0,p) 
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which gives, for p ^ S(Lfe, f), 

\dlS{r,k,p)\ < max \d^pF{r,k, p)\. 

tx Kj 

(Here we used that \dpUj{p)\ < x + ^- ) The higher order derivatives are estimated 
in the same way and this gives 

\dlS(r,k,p)\ < Ct.i(iV^^^)l^l maxja^F(r,fc,p)| 

for any |j| < s,, where Ct. is an absolute constant. 

By Lemma IA.2[ there exists a C'^-function : 2? —>■ K, being = 1 outside 
S(Lfc, k) and = 0 on S(Lfe, f) and such that for all j > 0 

\9k\a(v) < + 

K 

Multiplying S(r, k, p) with gk{p) gives a C®*-extension of S{r, k, p) from 2?\E(Lfe, k) 
to V satisfying the same bound as S(r, k, p). 

It follows now, by a classical argument, that the formal solution converge and 
that \dj,S{r, 9, p)\ and 0, p)\ fulfills the estimates of the lemma. When sum¬ 
ming up the series for 0 , p)| we get a term but the factor ^ 

disappears by replacing N by CN. 

By construction S and R solve equation (15.411 for any p G 'Di- □ 


5.3. The second equation. Concerning the second component of the homological 
equation we have 


Lemma 5.3. There exists an absolut constant C such that if (lot holds, then, for 
any N > 1 and 

0 < K < c', 

there exists a closed set (D 2 = T> 2 (/i, k, N) C T), satisfying 

Leh{V\V2) < 

Oo Oo 

and there exist -functions and Ry, on x 2? —>• Y^, real holomorphic in 9, 
such that for p € 2?2 

(5.7) du{p)Swi9,p) - A{p)JSy,{9,p) = fw{9,p) - Rw(9,p) 

and for all {9, p) G T^, x T), a' < a, and |j| < s* 


(5.8) 

(5.9) 


msy,{9,p)\\^ <c 


1 + 


mRy,{9,p)\\-, <C- 


„-{cr-cr’)N 




(tr — cr')” b.xr.D 

The exponent exp only depends on ^ and n = ffA, and C is an absolut constant 
that depends on c and sup^, |a;|. 


Proof. Let us re-write (15.71) in the complex variables ^ and 77 described in section 
13.21 The quadratic form {1 /2){w, A{p)w) gets transformed, by re = Uz, to 

{^,Qip)v) + ^{zj^,H'{p)zj^), 
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where Q' is a Hermitian matrix and H' is a real symmetric matrix. Then we 
make in (15.71) the substitution S = HJSw, R = *URw and F = Hjfw, where S = 
, S^), etc. In this notation eq. (15.71) decouples into the equations 

+ iQS^ = 

= F^ -R^ 

- HJS^ = F^ - R^. 


Let us consider the first equation. Written in the Fourier components it becomes 

(5.10) ((fc,w(p))/ + Q)5«(fc) = -i(F«(fc)-i?«(fc)). 

This equation decomposes into its “components” over the blocks [a] = [a] a and 
takes the form 

(5.11) Lk^[a]ip)Sla]ik) =: ((fc,w(p)) + Qla])Sla]ik) = -i(F[o](fc) - Rla]ik)) 

- the matrix being the restriction of to [a] x [a], the vector being is the 
restriction of F^ to [a] etc. 

Equation (15.101) has the (formal) solution 


Sla]ik,p) = 


-(Rk,[a](p)) ^iR[a](k,p) if |fc| < iV 


if not 


and 


For fc ^ 0, by Lemma 1412 


B (k o)- I 1^1 > ^ 

0 if not. 


\\iL,,[a]{p))-^\\<- 


for all p outside some set E(Lfejo], k) such that 


and 


^(F \ E(Lfcj„], k), E(Lfc.[a], |)) > 


'R2='R\ [J L;fcJa](K), 

o<|fc|<Ar 


fulfills the required estimate. For fe = 0, it follows bv 16.41 and 11.51 that 

||(Lfe.H(p))-'||<^<-. 

C t\ 

We then get, as in the proof of Lemma lOl that F[o](fc, •) and R[a]{k,-) have 
C'’*-extension to V satisfying 

P^FH(fc,p)|| < Ct.i(7V^)'^'' max 


K ' n 


and 


P^i?[a](fc,p)||<Ct.||a^F[,](fc,p)||, 
and satisfying (15.111) for p G 1^2■ 

These estimates imply that 

\\dj,S^{k,p)\\j < Ct.i(iV^^)'^' max ||a^F«(/c,p)||-y 

ri r\y 
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and 

Summing up the Fourier series, as in Lemma 15.21 we get 


P^^«(0,p)|l7<Ct.- 


;(cr — a')' 


K 0<Z<j |cje|<o. 


and 

\WpR^{S,p)\\^ <Ci.^ - - sup || 9 ^F’^(-,p )||7 

(cr-o-')" |ae|<^ 

for (0,p) G T^, X 2i>, 0 < cr' < ( 7 , and |j| < s*. This implies the estimates (I5.8p and 
dSSl) - the factor ^ disappears by replacing N by Ct.A^. 

The other two equations are treated in exactly the same way. 

□ 


5.4. The third equation. Concerning the third component of the homological 
equation, (lOl) . we have the following result, where for a solution Sww{(^,p) we 
estimate separately its mean-value SwwiO,p) and the deviation from the mean- 
value SyjyjiO^P^ P^ . 

Lemma 5.4. There exists an absolut constant C such that if (lea holds, then, for 
any iV > 1, A' > A > 1, and 

1 , 

«<pC, 

there exist subsets 'D^ = 'D^(h, k, N, A') C T>, satisfying 
Leb(P\p 3 ) < 

^0 ^0 

and there exist real -functions B^w '■ B fl AfJ^A' and Sww, Rww '■ x 

T) —>■ real holomorphic in 9, such that for all p G (D 3 

(5.12) d^(^p)Su,w{9, p) - A{p)JSww[9, p) + Sww{9, p)JA{p) = 

fww {9 , P^ Byjyj (p) Rww {9 , p) 
and for all {9,p) G T^, xT), a' < a, and |j| < s* 


(5.13) 


\diSww{e,p)\\, 




c(cr — <j') 






(5.14) ||d^i?„^(0,p)||^,^^<CA'A-P= 


„-{a-a')N _|_ g-(7-7')A' 

(a — cr')" 


fT 

\f , 

7,x,X> 


(5.15) 


KB,, 


.(p)|| <CA'A-P=|/^|.,. 


for any 7* < 7^ < 7. 

The exponent exp only depends on ^ and n = ffA. The exponent exp 2 only de¬ 
pends on d*, m*. The exponent a is a positive constant only depending on Sr, 

C is an absolut constant that depends on c and sup^j |a;|. 
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-’ww^ Rww 


and Bww verifying the 


Proof. It is also enough to find complex solutions Su 
estimates, because then their real parts will do the job. 

As in the previous section, and using the same notation, we re-write (15.121) in 
complex variables. So we introduce S = *US(^_qU, R = *URc^_qU, B = HjBi^^^U and 
F = In appropriate notations (15.121) decouples into the equations 

-b iQS^^ + - i?«, 

= F^'^ - B^^ - 




iQS^^ 


JH = - R^^ 


and equations for 5'’*’*, 5’’^^, Since those latter equations are of the 

same type as the first four, we shall concentrate on these first. 

First equation. Written in the Fourier components it becomes 

(5.16) {{k,io{p))I + Q)S^^{k) + S^^ikfQ = -i(i^««(fc) - 4,oB - R^^{k)). 

This equation decomposes into its “components” over the blocks [a] x [ 6 ], [a] = [a] a, 
and takes the form 


(5.17) L(fc, [a], [&],p)^[^](fc) =: (fc,w(p)) 5[^](fc) -b (5[a](p)^{^](A:)+ 


S^]ik) ^Qwip) = -KFi:i{k,p) - i?p(fc) - 4,oi3[:j) 

- the matrix Q^a] being the restriction of to [a] x [a], the vector being the 
restriction of F^^ to [a] x [5] etc. 

Equation (15.161) has the (formal) solution 

eWifc oj = / if W) < and |fc| < N 

[“ 1 ^ 1 0 if not , 


■fb] 






0 


Kik.p) = 


For fc 7 b 0, by Lemma l4A 


_ / Faik, p) if dist([a], [^]) > A' or |fc| > N 


if not. 


\\{Fk.[a],[b](P)) ^11 


for all p outside some set E^, [„]_[{,] (k) such that 


and 


dist(2^\ Sfe,[a],[h](K),Efej„]j,,](-)) > 

F3=F\ [J Sfc_[a]Jb](K) 

0<|fe|<Af 

[a], [6] 

fulfills the required estimate. For fc = 0, it follows bv 16.41 and 11.51 that 









32 


L. HAKAN ELIASSON, BENOIT GREBERT, AND SERGEI B. KUKSIN 


We then get, as in the proof of Lemma [Ol that •) and •) have 

C®*-extension to V satisfying 

and 

p^<(fc,p)||<ct.||a^A(fc,p)||, 

and satisfying (15.171) for p £ 1 ) 3 . 

These estimates imply that, for any 7 * < 7 ^ < 7 , 

A exp„27(iA V -I- A I -I 

\\dj,S^Hk,p)\\ 6 iY,„Y^,)<Ct.A' --- {N^y^\maxMdlF^Hk,p)\\BiY^,,Y^,) 

tv tv Ut7 

and 

\\dj,R^Hk,p)\\BiY^„Y,,) < Ct.A'A^-^\\d^^F^Hk,p)\\BiY^„Y^,). 

The factor occurs because the diameter of the blocks < cIa interferes 

with the exponential decay and influences the equivalence between the Z^-norm and 
the operator-norm. The factor A'A®’'p occurs because the truncation < A' -|- cZa of 
diagonal influences the equivalence between the sup-norm and the operator-norm. 
These estimates gives estimates for the matrix norms and, for any 7 * < 7 ' < 7 , 

Aexp 27(iA V -I- ri I - I 

||a^5«(fc,p)|k.<Ct.---max 

rC rv 

and 

\\dj,R^Hk,p)\U,.<Ct.\\df,F^Hk,p)\U,.. 

Summing up the Fourier series, as in Lemma 15.31 we get that S^^{d,p) satisfies 
the estimate (I5.13|) and that R^^{9,p) satisfies the estimate (I5.14|) . 

The third equation. We write the equation in Fourier components and decom¬ 
pose it into its “components” on each product block [a] x [Zi], [ 6 ] = F-. 

L{k, [a],[Z)],p),s|^j(fc) := {k,u}{p)) (fc)-h Q[a](*)- 

i5£j(fc)JiL(p) = -i(Ff](*>P) - - Rnik)) 

- here we have suppressed the upper index ^zj 7 . 

The formal solution is the same as in the previous case and it converges to 
functions verifying (15.131) . (I5.14|) and (I5.15F by Lemma 031 and by (II.6F 
The fourth equation. We write the equation in Fourier components: 


L{k, [a], [6],p)5£j(fc) := (fc,u;(p)) - iH Jip)S^^ik)+ 

(A) jij(p)=-i(^':j(fc,p) - 4.0 a'3 - 

where [a] = [6] = - here we have suppressed the upper index zj^zj^. 

The equation is solved (formally) by 




-L(fc,[a],[Z,],p)-iiF[5'(fc,p) if 0 < |fc| < A 

0 if not. 


Rt!]^k,p) 


F[g(fc,p) if|fc|>iV 
0 if not; 
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and 

The formal solution now converges a solution verifying (I5.13L (15.141) and (15.151) 
by Lemma l4^ 

The second equation. We write the equation in Fourier components and decom¬ 
pose it into its “components” on each product block [a] x [5]: 

L(fc, [a], =: {k,w{p)) -b Q[a](p)*5[^](fc)- 

s^]{k)Q[t]{p) = -KF^iKp) - - 4 . 050 ) 

- here we have suppressed the upper index ^rj. This equation is now solved (for¬ 
mally) by 

with 


o') = J 5(fc, [a], [b],p) 4F0(A:,p) if dist(H, [b]) < A' and 0 < |fc| < iV 
[“1^ 1 0 if not, 


?bfu / K{k,p) if dist(H, M) > A' or |fc| > N 


Kik,p) = 


if not. 


and 


= / ^a{0,P) if dist([a], [6]) < A'and fc = 0 
10 if not. 


We have to distinguish two cases, depending on when fc = 0 or not. 
The case k ^ 0. 

We have, by Lemma 1431 


ll(5fe.H.M(p))-'ll 
for all p outside some set Efc_[a],[b](«^) such that 


1 

< - 

K 


dist((D \ Efe [Q]_[{,](«), Sfejaj.fb] ( 2 )) ^ ct. i 


and 


F>3=V\ [J Efcja],[(,](K) 

0<|fe|<Af 

[a].[b] 

fulfills the required estimate. 

The case k = 0. In this case we consider the block decomposition and we 
distinguish whether |a| = |6| or not. 

If |a| > |6|, we use (16.41) and (11.61) to get 


\a{p) - /3(p)l > c' - 


6 S ^c' ^ _ 
(a)^ ~ 1^ 


This estimate allows us to solve the equation by choosing 
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and 

with 




-1 p[b] / 


i'W/'n ^Mi ^ /xrKdL^'ib'l pii 


||a^^f:j(0,p)|| < Ct.-(7V:^)l^l max 
which implies (I5.13|) . 

If |a| = |6|, we cannot control \a{p) — /3(p)| from below, so then we define 


5£](0) = 0 

and 

i?^(p) = F,'(0,p)), i?^(0)=0 for Ha'= [6]a' 

Rl{Q,p) = Ft{0,p) Bl = Q, for [a]A'[6]A'. 

Clearly R and B verifiy the estimates (I5.14|) and (15.151) . 

Hence, the formal solution converges to functions verifying (j5.13ll . (15.141) and 
(15.151) by Lemma [4.31 Moreover, for p G V', these functions are a solution of the 
fourth equation. 

□ 


5.5. The homological equation. 

Lemma 5.5. There exists a constant C such that if dsa) holds, then, for any 
N > 1, A' > A > 1 and 

1 , 

»<pC, 

there exist subsets T)' = k,N) C T>, satisfying 

LehiV\V') < 

Oo ^^0 

and there exist real jet-functions S^RG f^) ct'^d /i+ verifying, for p G T>', 

(5.18) {h,S} + p = h+ + R, 

and such that 

h + h+GMR^{A’,5+) 

and, for all 0 < a' < a, 


(5.19) 

(5.20) 

(5.21) 




'f ,>c,T> 




<ix(iv 2 ^)-|/^'|., 




'y,>c,T> 




7 


; 

7,>f,X> 


for 7 * < 7^ < 7 , where 

X = max(l,^2)^ 
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The exponent expj^ only depends on ^ and The exponent exp 2 only de¬ 

pends on d^:,m^ and 4j^A. The exponent a is a positive constant only depending on 
C is an absolut constant that depends on c and supp |a;|. 


Remark 5.6. The estimates (15.191) provides an estimate of i5-|-. Indeed, for any 
a,6 G [a]A' 


< Ct.(A'r \r \ 


C 

Since #[a]A' < (A')®’'p we get 

||5^i3(p)[a],,||<Ct.(A'rP|/^|..M 
This gives the estimate of <5+ — <5. 


(a)^' 


7 ,’x,r> (a)^' 


Proof. The set V will now be given by the intersection of the sets in the three 
previous lemma of this section. We set 

h+{r,w) = fr{r,0) + ^{w,Bw) 


S{r,6,w) = Sr{0,r) + {Snj{0)w) + ^{Si^n,{d)w,w) 

and 

R{r, 9, w) = Rr{r, 9) + {R^{9),w) + ^{Rww{d)w, w). 

These functions also depend on p gV and they verify equation (15.181) for p gV . 
If X = (r, 9, w) G (a, p), then 



|h+(x)| < 

1 fTi 

J K.A* 

7,>f ."D 


Since 






— 1 Z3(T7. .'Ey. ) 

it follows that 





|h+(x)| 

-G 

o 

VI 

1 o-.Ai max(l,/i^). 

j.x.T) 


We also have for any x = (r, 9, w) G 0-y'{a, p), "i* < 7, 


II Jdh-|_(x)||7' < Ct. 1/^1 <y,n + ||i?w||7'. 

Since 

||B||,,^>||i?||,,^>||S|lB(V,,.V,,) 

it follows that 

II Jdh+(x)||.^/ < Ct. I/'^I cr./i max(l,^). 

'y,>c,T) 

Finally Jd^h+{x) equals JB which satisfies the required bound. 

The estimates of the derivatives with respect to p are the same and obtained in 
the same way. 

The functions S{9, r, () and R{9, r, f) are estimated in the same way. 
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5.6. The non-linear homological equation. The equation (15.2p can now be 
solved easily. 

Proposition 5.7. There exists a constant C such that for any 

S<^c\ 

and for any 

N>1, A' > A > 1, K<^c' 

o 

there exists a subset T)' = 'D{h, k, N) C T>, satisfying 
Leh{V\V') < 

Oq Oo 


and, for any f € g,,!)), fi <1, 

s = If’^l <y,u and f, = \f\o-,u , 

there exist real jet-functions S, R € g,) and h+ verifying, for p € V, 

(5.22) {h,S} + {f-f,S}^ + f = h+ + R 


and such that 

h + h+ €j\fR,,{A',S+) 
and, for all a' < a and p,' < p, 


(5.23) 


\h+l,^^ <CXYe 

7,>f,X> 


(5.24) 


7,>f,T) 


< C-XYe 

K 


(5.25) \R\Y,t, < C ATe, 

7 


for 7 * < 7 ' < 7 , where 


and 


^ NA’ey^- 
Y = (- X + 


The exponent exp^ only depends on ^ and ffA. The exponent exp 2 only depends 
on d^:,mt:,St, and ffA. The exponent a is a positive constant only depending on 
C is an absolut constant that depends on c and supp |a;|. 


Remark 5.8. Notice that the “loss” of S with respect to k is of “order” 4s* -|- 3. 

However, if y, 6 and £, = \ f\ (y,v are of size < k, then the loss is only of “order” 1. 

Y ,k,t> 
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Proof. Let S' = + 5'i + S '2 be a jet-function such that Si starts with terms of 

degree f in r, w and S 2 starts with terms of degree 2 in r,w - jet functions are 
polynomials in r, w and we give (as is usual) w degree 1 and r degree 2 . 

Let now cr' = CT 5 < 174 < (73 < (T 2 < <ti < co = o' be a (finite) arithmetic 

progression, i.e. Uj — cFj+i do not depend on j, and let and ^' = /is < ^4 < /is < 

M 2 < Ml < Mo = M be another arithmetic progressions. 

Then {M', S} -I- {/ — /^, S}"^ -\- f'^ = hj^-\- R decomposes into three homological 
equations 

{h',So} + f = {h+)o + Ra, 

{h', Si} + fl = (/i+)i + i?i, /i = {/ - f, So}, 

{h', S 2 } + /J = (h +)2 + S 2 , /2 = {/ - f: Si}. 

By Lemma 15.51 we have for the first equation 

|(h_|_)o|cri,Ai ^ JCe, |i?o|o'l,M <XZ£, 

1 ' 


[Solcu.M < —XYe 

Kj 

where 

— a' 

and where Y, Z are defined by the right hand sides in the estimates (15.201) and 

(ICT) . 

By Proposition 12.51 we have 


Cl = 1 / 1 1-2,< -XYWfe 

7,x,r> K 


where 


W = C{ 


(a-a') (m-mO 


n)' 


By Proposition [ 73 ] £1 = II- 2 . 1*2 satisfies the same bound as Ci 




By Lemma 15.51 we have for the second equation 


KMil -3.A‘2 ^ ^£ 1 , 

'y.>c,T> 


■y' ,>c,T) 


^ —XXsi. 

i,>c,v n 


By Propositions 12.41 and 12.51 we have 


6 = I/2I-4.M4 < -XFfPCiei, 


and £2 = I/; 


- 4 ,m 4 satisfies the same bound. 

7,>f.'D 

By Lemma 15.51 we have for the third equation 


I (h_|_)2 |-5,M4 < ^£2, 

'y.^c.'D 


|i?2|-5,/44 <XZ£2, 

-y' ,>c,T> 


15^2 — XXS2- 

l,>c,V K 

Putting this together we find that 

e + ei + £2 < (1 + -XYWffe = Te 
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and 

<XTe, <XZTe, 

7',>c.D 

< -XYTe. 

■rX.'D « 

Renaming X and Y gives now the estimates. □ 

6 . Proof of the KAM Theorem 

Theorem 13.61 is proved by an infinite sequence of change of variables typical for 
KAM-theory. The change of variables will be done by the classical Lie transform 
method which is based on a well-known relation between composition of a function 
with a Hamiltonian flow and Poisson brackets: 

from which we derive 

fo<^s = f + {LS} + [\i-s},s} o<i>*s dt. 

Jo 

Given now three functions h, k and /. Then 
(h + k + f) o <I)g = 

h + k + f + {h + k + f,S}+ [ (1- t){{h + k + f, 5}, S} o dt. 

Jo 

If now S' is a solution of the equation 

(6.1) {h,S} + {f-f,S}^ + f = h++R, 
then 

{h Y k Y f) o = h Y k Y Y /+ 

with 

(6.2) U = RY{f-nY{kY f, S}Y{f- f, S}-{f- f, 

+ [ {l-t){{hYkYf,S},S}o<^ldt 

Jo 

and 

(6.3) fl = R+{kY S}^ + {[\l- t){{h YkYf, S}, S} o dtf. 

Jo 

If we assume that S is “small as” /^, then /J is is “small as” fc/^ - this is the 
basis of a linear iteration scheme with (formally) linear convergence. E3 But if also 
k is of the size f'^, then /+ is is “small as” the square of - this is the basis 
of a quadratic iteration scheme with (formally) quadratic convergence. We shall 
combine both of them. 

First we shall give a rigorous version of the change of variables described above. 


it was first used by Poincare, credited by him to the astronomer Delauney, and it has been 
used many times since then in different contexts. 
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6.1. The basic step. Let h G S) and assume x > 0 and 


(6.4) 


A < 

- C ’ 


where C is to be determined. 

Let 

7 = (7,m*) > 7* = (0,m*) 
and recall Remark 1 5. II Let A^>1, A'>A>1 and 

1 , 

K < -C . 

Proposition 15.71 then gives, for any / S 71y,>c,z)(o',/r), /r < 1, 

£ = \ f'^\ and £ = |/| o-,/^ , 

a set V = V'{h,K,N) C V and functions S,h+,R satisfying (I5.23|) + ()5.24|) -I- (I5.25I) 
and solving the equation (EH), 

{h, S} + {f- f, 5}^ + f = h+ + R, 

for any p G V. Let now 0 < a' = a 4 < as < a 2 < ci < ao = a and 0 < p' = pi < 
P 3 < P 2 < Pi < Po = P (finite) arithmetic progressions. 

The flow 4>g. 

We have, by (15.241) . 

IS'lo-ij/ii <Ct.—XY£ 

7 ,xr,X> K 

where X,Y and Ct. are given in Proposition 15.71 i.e. 

^ = (t-V- " = (t -A?-n) ^ 

[as - ai){po - pi) [a - a'){p- p') 

Y ^ j- X + (5 + 

K 

- we can assume without restriction that exp 2 > 1 . 

If 

(6.5) . 4 ^, 

and C is sufficiently large, then we can apply Proposition l2.7l i') . By this proposition 
it follows that for any 0 < t < 1 the Hamiltonian flow map 4>g is a C®*-map 

0^'{ai+i,Pi+i) xV ^ 0^'{ai,pi), V 7 * < 7 '< 7 , i = 1,2,3, 

real holomorphic and symplectic for any fixed p G T>. Moreover, 

|l4($4x,-)-4ll7' <ct.iAr£ 

rv 


and 


||4(d$4(r,.)-/)|L ,^<Ct.iAF£ 

^ ' rxj 


for any x € 0^'{a2, P 2 ), 7 * < 7 ^ < 7 , and 0 < |j| < s*. 

A transformation. 

Let now k G h) and set 

p = |/c| cr,fj, . 
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Then we have 

(/i + fc + /) o = h k /i_|_ + 
where /+ is defined by (j6.2p . i.e. 

U = R + {f-f) + {k + f, s} + {f- f, s}-{f- f, s}^+ 

+ [ {l-t){{h + k + f,S},S}o<^*gdt. 

Jo 

The integral term is the sum 

[\l-t){h+ + R-f,S}o<i>%dt+ [\l-t){{k + f,S}-{f-f,S}^,S}o<^%dt. 
Jo Jo 

The estimates of {k + /^, S} and {/ — /^, S}. 

By Proposition I2.5f i') 


Hence 

( 6 . 6 ) 

Similarly, 

(6.7) 


{k + f'^, S}\<72,fJ.2 < Ct.^ |S'|cri,/xi Ifc + /^|cri,/xi 

•y,>c,T> 7,Q:,X> 

|{fe + /^,S'}U2.M2 < Ct.-X^Y{r] + s)e. 

'y,a,T> K 

\{f - f, S}\<T2,n2 <Ct.-X^Y{f + e)e. 

'y,a.,T> K 


The estimate o/ {/i+ — S'} o <1)^. 

The estimate of /i+ is given by (15.231) : 

|/i+|o-i,Aii < Ct.XYe. 

This gives, again by Proposition 12.Sl iil. 

\{h+- f^,S}\a2^n2 <Ct.-X^Y^e^. 

7,Q:,'D K, 

Let now F = — /^, Sj. If e verifies (16.51) for a sufficiently large constant C, 

then we can apply Proposition 12.Tf iil. By this proposition, for |t| < 1, the function 
F o e T},>r,-D(o- 3 , Ms) and 

( 6 . 8 ) \{h+- f,S}o^Ua2.n3 <Ct.-X^Y^e^. 


The estimate of {R, S} o <I>g. 

The estimate of R is given by (15.251) : 


\R\<T 2 ,ni < Gt.XYZ^e, 

7^ 

where 

Zy = . 

Then, as in the previous case, 

(6.9) |{S,S}o$‘g|<T3,/.3 < Ct.-X3r2v £2. 

The estimate of {{/c + /, Sj — {/ — /^, S}^, Sj o 
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This function is estimated as above. If T" = {{k +f, S} — {f — f'^, S}'^, S}, then, 
by Proposition 12.41 and Proposition 12.5l ii. 

<Ct.{-X^Yy{r] + Oe^ 


7,q;,X? 


and by Proposition 12.7f i 


,1 


(6.10) {{fc + f,S}-{f- }^ 5} o a 4./.4 < Ct.i-X^YYir, + 

K, 

Renaming now X and Y and replacing N by 2N now gives the following lemma. 

Lemma 6.1. There exists an absolute constant C such that, for any 

1 


h & J\fX x> Q, 5 <—c', 

o 


and for any 


N >1, A' > A > 1, K< —c', 

o 

there exists a subset V' = ^{h, k, N) C T), satisfying 
Leh{V\V') < 


and, for any f £ /r), yi<l, 

(T 


e= / a.M and ^ 

7,x,P 


satisfying 


£ < 


1 K 


Y = ( 2 dA±i)exp 3 ^ 


and for any k £ /r), there exists a C®* mapping 

f 

,k-- n' )’ V74,<7'<7, 


$ : e>y(CT',/r') xV-^Oyia- 

ilomorphic and symp 
f+,R+ £ ,n') and 


a — o' T~ h-' 


^7 ^ J ^ ^ ^ 2 2 

real holomorphic and symplectic for each fixed parameter p ^ T>, and functions 


such that 
and 


h + hj^ € J\f J' 

(^h-\-k-\-f)o^ — h-\-k-\- hj^ + /-I- + -R+, yp G 


\h+\Yy +|/+-/L'y <CAre, 

7 .>c.X> 'y.>c,'D 

/II.',/ <C-XYi\k\a,^ + ^e-^-Y)N ^ 


7,>f,X> 


^.>c.T> 


and 


<CXYe-^^-^'^^'e 

7^>f,X> 


for any 7 * < 7 ^ < 7 . 
Moreover, 


1 


\m<^{x,p) - ai)||y + \md<^{x,p) - /) I , ^ < C-XYe 


~ K 


for any x £ , p'), 7 * < 7 ' < 7 and \j\ < s*, and for any p £7). 
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Finally, if p = (0, p 2 , ■ • ■, Pp) and p) = 0 for all p, then /+ — / = R+ = 
/i+ = 0 and $(a;, ■) = x for all p. 

The exponent exp^ only depends on ^ and ffA. The exponent exp 2 only depends 
on d,t,m^:,s^, and ffA. The exponent expg only depends on s*. The exponent a 
is a positive constant only depending on C is an absolut constant that 

depends on c and sup-p |a;|. 

6.2. A finite induction. We shall first make a finite iteration without changing 
the normal form in order to decrease strongly the size of the perturbation. We shall 
restrict ourselves to the case when N = A'. 


and for any 


Lemma 6.2. There exists a constant C such that, for any 
hGAfF^{A,S), X > 0, 5 <—c', 

o 

A' > A > 1 , K < 

o 

there exists a subset V = ^{h, k, A') C V, satisfying 

Leh{V \V')< C(Ay^P^ (1^) 

Oo 00 


and, for any f £ /r), p<l, 

(T 


satisfying 


£ < 


e = / a.M and ^ 


1 ^ I ^ = ((^^^5^1ogirP=, a'<a,p'<p 
■ 


C XY" [ Y = ( X+^+^ )exp3^ 
there exists a C'** mapping 

$ : Oy{a',p) xV ^ Oy(a- ^ ^ Vy* < 7' < 7, 

real holomorphic and sympleetic for each fixed parameter p £ T), and functions 
f £ Ty,^,v{A,p') and 

h' £ MX^{A', 5 '), 

such that 

(/i + /)o$ = /!' + /', VpeP', 

and 

\h'-h\y^y <CXYe, 

7 ,>tr,'D 

e = \fV,F <^ + CXY£ 

and 

e' = \{ff\a'X < CXY{e-i^‘^-'^'^^' + e- 5 (T'-T'')^')e, 

7^,x,X? 

for any 7* < 7' < 7. 

Moreover, 

||a^($(x,p) - x)||y + \\dl{d^{x,p) - /)|| ^ < C-XYe 

for any x £ Oy'{a', p’), 7* < 7' < 7 and \j\ < s*, and for any p £ 7 ). 
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Finally, if p = (0, p 2 , ■ ■ ■, Pp) and p) = 0 for all p, then f' — f = h' = 0 and 
$(x, ■) = X for all p. 

The exponent expj^ only depends on ^ and ffA. The exponent exp 2 only depends 
on d,t,m^:,s^, and ffA. The exponent expg only depends on s*. The exponent a 
is a positive constant only depending on s^,^,^- C is an absolut constant that 
depends on c and sup^j |a;|. 


Proof. Let N = A'. Let cri = cr — , pi = p — and ax+i = A , pk+i = 

and let and be arithmetical progressions. We take K such that 


i.e. K < ((j-(T')A'(logf)-b 

We let /i = / and fei = 0, and we let ei = [f^] ^ = £, = [fi] ^ 

(5i = i5 and pi = [fci] a.p = 0. 

Define now 




7,>f,'D 


with 


£j-\-l — C — XjYj[Tlj-\-£i-\-6j)£jj 

Ki 

^ 4+1 = Vj+i = Vj + C^AjYjEj, 


X- = ( _ AA e'^ _ _ , x + ^ + f,j xexp. 


where C, exp 2 , expg are given in Lemma l6.ll One verifies by an immediate induction 
that 


Sublemma. There exists an absolute constant C such that if 

^1 K 

then, for all j > 1, 

1 ^ x'^y'^ 

£j < ^^2^ and Ej < (Ct.^^ei)^-i£i, 

{fj - 6) + iVj - m) < C'XiYiEi. 

The constant Ct. only depends on C and expg. 

We can then apply Lemma 16.11 K times to get 

: Oy{aj+uPj+i) xV ^ Oyiaj - 7 * < 7' < Ij 

and and Rj+i such that, for p G V, 

(h + kj + fj + Sj) o = h Y kj Y + /j+i + Aj-i-i + Sj o <1?^ 

with kj+i = kj Y hj+i, ki = 0, = Rj+i Y Sj o <1)^, Si = 0. 

We then take $ = ihi o • • • o <1)^, h' = hY k^+i and /' = fx+i + >5'a+i- 
Then 

YC'XiYiE, S' <C{A'f^P^XiYiE 

7,>f,X’ 


e = \f\Y.p' <f + C'XiYiE 

^ .>c.T> 


and 
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and, for p gV', 

/ ^ tv 

For the estimates of write o • • • o and = id. For (x, p) € 

xV we then have 

K 

||$(x,p) - xIIy < W'^jix.p) - ^'j + i(x,p)||y. 

4 = 1 


Then 


||«'j(x,p) - 'iij+i{x,p)\\Y = \\<i>j{'i/j+i{x,p),p) - 'iij+i{x,p)\\Y 
is 

< Ct.-XiYiSj max(-----), 

K k-cr'l |/i-/i'| 

by a Cauchy estimate. Hence 

||$(a:, •) - a;||y < Ct.-XiFi max(-—77)e. 

K k-cr| Im-mI 

The derivatives with respect to p are obtained in the same way, as is also the 
estimates of d^. 

The result now follows if we take C sufficiently large and increases the exponent 
exp 2 . 

□ 


Proof of sublemma. The estimates are true for j = 1 so we proceed by induction 
on j. Let us assume the estimates hold up to j. Then, for k < j, 


+ ^ + 6 + ^C'XiYiEi 
“ K 




and 

£j+i < 2®’'P3:^iF[2C'Xiyiei +ei +£i\ej < Ct.—!-^£i£j. 

tv tv 

Then 

0+1 ^ 0 + ct.XiYi(£i + • • • + £j+i) < 0 + 2Ct.Xiyi£i 
and similarly for Pj+i- 


6.3. The infinite induction. We are now in position to prove our main result. 
Theorem 13.61 

Let h be a normal form Hamiltonian in and let / € Ty^j.c,T>io', p) be 

a perturbation such that 


0<£ = f 


■y.x.D 


We construct the transformation <i> as the composition of infinitely many transfor¬ 
mations as in Lemma 16.21 We first specify the choice of all the parameters for 
J > 1- 

Let C, expj, exp 2 , expg and a be the constants given in Lemma 16.21 
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6.3.1. Choice of parameters. We assume (to simplify) 1 and A > 1. By 

decreasing 7 or increasing A we can also assume 7 = (cIa)”^. 

We choose for j > 1 

A j /I 15 

M4 = (2 + ^)/^ and aj = {- + —}a. 

We define inductively the sequences Sj, A^-, Sj and £,j by 


( 6 . 11 ) 


where 


Bj+i = e^^CXiYie ei = e 

Aj+i = AKj max(^- 3 ^, dA,) log ^ Ai = A 

7 j+i = (dA,+i)"^ 71=7 

Sj+i = Sj + CXjYjSj Si = S>0 

. ^i+i “ ^ A 




W = (^^- 1 -r log i)®='P 2 = ( 

Y = / XT’S?+^ 7 'jexp, 

3 \ Ki ' 


A,-A,-+ie4^+^ 

(Tfl 


log i)®XP2 


The Kj is defined implicitly by 


1 Ki 


S') —-, 


These sequences depend on the choice of Kj. We shall let Kj increase like 

Kj = K^ 

for some K sufficiently large. 

Lemma 6.3. There exist constants C and exp' such that, if 

K> C 


and 


then 

(i) 


(ii) 


F < 1 ( O'M \exp'/ d N,exp3 , 

~ C'^KAlog^ + d + 

00 . 

J^CXkYkek < 2CXiYie < ^c'. 

k^l 


00 

C max( 

fc=i 


1 1 

Cfc — CTfc+i ’ fj.k — fJ-k+1 


)XkYkek < 5C'max(-, -)Aiyie < 

(7 


1 


/ 


C . 


(iii) 


^CA-7(i 

1>1 


-l!l^i)(^A)“ < ^exp' 

So ^ So ~ ap 


oo So 


C is an absolut constant that only depends on fS, x, c and supp |a;|. The exponent 
exp' is an absolute constant that only depends on /3 and k. The exponents a' and 
/? are positive constants only depending on St,, 
















46 


L. HAKAN ELIASSON, BENOiT GREBERT, AND SERGEI B. KUKSIN 


Notice that (i) implies that 


Proof. Aj+i is equal to 

4i4rj max(- - -,dA, )log- < (Ct.-log i)(2N:)'’^ A“ = A(2Ky^ A'^, 

(7j (Tji_j_2 S (7 S 

which, by an induction, is seen to be, by assumption on e, 

< M(27^)“A)“' < (i)“' 

£ 

if a is, say, at least 6. In the same way one sees that 

Xj < (i) 2 exP 2 a\ 

(i). For j = 1, (i) holds by assumption. Indeed, by definition 
(CXiyiei)i+'^"P3 = «;1+-P3 = CAiFi^fP^ei < h ‘^"P^(x + J + 0^"P3£, 

^ C7/i ^ 

which is 

- Uc ^ 

by assumption on e. 

Assume now (i) holds up to j — 1 > 1. Then 5j < 5 + 2CXxYie and fj < f + 
2CXiYie^ and hence 

F, < (X±i±i±i^iyi^)exp3 < ct.yi(^rP3, 

Kj Kj 

and, by the definition of Kj, 

^i+exp3 ^ cXjYjejKy""^^ < Ct.FiKf < XjS^^ 

by assumption on e. Hence 

CXjYjSj = Kj < < e2bif.-2exp3 A < bKj ^ ^ — I —^ 

3 - - ’ 2 (exp 3 + 1 ) ’ 

if K is large enough - notice that j > 2. This implies that 

j 

CX^YhSk < 2£^^3 < £ < C'AiFi£i 

k=2 


if K is large enough. 

The proof of (ii) is similar. To see (hi) we have for j >2 


A^^'Pi K-“ 

^j+i S' 


(A®(^P'/tj)“ < (Ayp'«;j)“ < (e2&R^-2(exp'+l)exp3a 




CX. 


which is 


< £^^J“ 


if K is large enough. 

Therefore 

Y < 2 A“Pik“ 

i>i 

if K is large enough. 


□ 
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6.3.2. The iteration. 


Proposition 6.4. There exist positive constants C, a' and exp' such that, for any 
h S S) and for any f G p), 0 < 7 , tr, // < 1 , 

I sT I 
£ = / 


*/ 


and 






< 1 / max(7 \dA) x-exp'. c' 

®£ ap. X + '^ + C 


expj^, 


then there exist a set T)' = ^'{h, f) C "D, 


Leb(P \ A) < C{ log )«(£)“', 

eafi do do 

and a C®* mapping 

$ : O-y, ((t/ 2, /i/2) X 2? —>■ (a, fi), 

real holomorphic and symplectic for given parameter p £ T), and 

c' 


such that 
verifies 


h' &J\fF^{oo,d'), d'< 


{h + f)o^ = h' + f' 


2 ’ 


and, for p G V, (/')^ = 0. 
Moreover, 


and 


\r-f\aM2 <C' 

7* 


|/l'-/iU/2W2 <C' 

7* 


for any x G O(o,m,){o'', h-') o^nd |/| < s*, and for any p GV. 

Finally, if p = (0, P 2 t ■ ■, Pp) and f'^{-,p) = 0 for all p, then h' = h and $(x, •) = 
X for all p. 

C is an absolut constant that only depends on /3, k, c and sup^? |a;|. The exponent 
exp' is an absolute constant that only depends on (3 and >c. The exponent expg only 
depends on s*. The exponent a' is a positive constant only depending on 

Proof. Assume first that 7 = 

Choose the number pj, Uj, £j, Aj, 7 ^, dj, , Xj, Yj, kj as above in Lemma [H31 with 

K = C. By the assumption on e we can apply Lemma 16.31 

Let hi = h, fi = f and 2?i = T). Lemma 16.31 now implies that we can apply 
Lemma im] iteratively to get for all j > 1 a set Tdj j,.i C Vj such that 

Leb(I), \Pj+i) < 

On On 
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a C®* mapping 

^ V 7 * < 7 ' < 

real holomorphic and symplectic for each fixed parameter p, and functions /j+i € 
(^j + l;/^i+l) 

^ J\f J~ jtc ) 

such that 

{hj + fj) o = hj^i + fj+i, Vp G T^j+i, 

with 


and 


fl^i\ 


<^}+i 

13 + 1,^■.'D 


< £j+l 


Moreover, 

and 


l/j + l kj+i.Mi+i <0+1- 

73 + 1 ,>f,D 

7j+i,>c,D 


P'($,+i(x, •) - x)||y + \\dl{d<t>,+3{x,-) - 7)11 , ^ < C-X,Y,e, 


for any x G O7'(cTj+i,/ij+i), 7* < 7' < 7^+1 and |/| < s*. 

We let h' = lim/ij, /' = lim/j and d) = $20- • •o$30.... Then (/i+/)o<f> = h'+f 
and h' and /' verify the statement. The convergence of 3 > and its estimates follows 
by Cauchy estimates as in the proof of Lemma l 6 . 2 l 

The last statement is obvious. 

If 7 > (dA)~^, then we can just decrease 7 and we obtain the same result. If 
7 < then we increase A and we obtain the same result. □ 


Theorem 13.61 now follows from this proposition. 


7. Examples 

7.1. Beam equation with a convolutive potential. Consider the d* dimen¬ 
sional beam equation on the torus 

(7.1) Mtt-I-A^ m-I--I-£ p(a;, u) = 0, xGT*^*. 

Here p is a real analytic function on T'^* x I, where 7 is a neighborhood of the origin 
in K, and the convolution potential V : T'^* —?► K is supposed to be analytic with 
real Fourier coefficients E(a), a G Z'^*. 

Let A be any subset of cardinality n in . We set L = lA' \A,p= (14)aeAj 
and treat p as a parameter of the equation, 

P = (Pai, . . .,Pa„) G 77 = [Pa[,Pa';] X •• • X [pa^ , Pa" ] 

(all other Fourier coefficients are fixed). We denote Pa = -I- V{a), a G Z'^*, and 

assume that Pa > 0 for all a G M, i.e. |a|^ + p'^ > 0 if a G A. We also suppose that 

pi A dll A dl2 G A, l\ yf l2- 

Denote 

= {a G £ : Pa < 0}, PI =: A, Coo = C\fF, 


7i+i, 




KAM FOR THE NON-LINEAR BEAM EQUATION 2 


49 


consider the operator 

A=|A^-hy* 1^/^ = diag{Aa,a e = 

and the following operator linear over real numbers: 

r aG£oo, 

\ oGT", 



Introducing the complex variable 

^ = ^(Ai/ 2 y _ zA-i/^u) = (27r)-'^/2 ^ , 

aGZ-i. 


we get for it the equation (cf. [3 Section 1.2]) 

( 7 . 2 ) i = i(h*f + (^) ) - 

Writing ijja = {ua + iva)!'/^ we see that eq. is a Hamiltonian system with 

respect to the symplectic form ^ dva A dua and the Hamiltonian h = hup + eP, 
where 

P = y^ G A~^/^ dx, duG{x,u) = g{x,u), 

and hup is the quadratic Hamiltonian 

hup{u,v) = AallAaP + + XI ■ 

aeA \ P / \ F / 

Here uj- = *'{ua, a G P) and H is a. symmetric 2N x 2 A-matrix.The 2N eigenvalues 
of the Hamiltonian operator with the matrix H are the real numbers {±Aa, a € P}. 
So the linear system (EH) ls=o is stable if and only if n = 0. 

Let us hx any n vector / = {/„ > 0, a G Al} with positive components. The 
n-dimensional torus 

f |i/)up = Ja, a G Al 
[ V'a = 0, a G C = Z"^* \ Al, 

is invariant for the unperturbed linear equation; it is linearly stable if and only if 
N = 0. In the linear space span{'(/;a,a G Al} we introduce the action-angle vari¬ 
ables {ra,9a) through the relations ijja = \J{la + Ta)e®®“, a G Al. The unperturbed 
Hamiltonian becomes 

hup = const -f (r,w(p)) + ^ Aal^-aP , 

^ 2 \ / 


with L>j{p) = (uJa = Aa, a G A), and the perturbation becomes 


P = e 


G(x,u(r,0;()(x))dx, u(r, 9; ()(x) = A 


i.e. 


= E 


aeA 


\/(dg + Xg) (6^^“ (fig + e 


E 

ae£ 


'0a£a + '>Pg‘P-g 

V^Ag 


In the symplectic coordinates {{ug,Vg), a G C) the Hamiltonian hup has the form 
(inj, and we wish to apply to the Hamiltonian h = hup + eP Theorem 13.61 and 
Corollarv l3.7l The assumption Al with constants c, c' of order one, Pi = P 2 = Ps = 
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2 holds trivially. The assumption A2 also holds since for each case (i)-(iii) the second 
alternative with a;(p) = p is fulfilled for some So ^ 1. Finally, the assumptions R1 
and R2 with x = \ and suitable constants 71 , cr, /r > 0 and 72 = m, are valid in 
view of Lemma 3.2 in [5]. More exactly, the validity of the assumption R1 is a part 
of the lemma’s assertion. The lemma also states that the second differential JcPf 
defines holomorphic mappings 

J(ff:OY{a,p)^M°,, 7 '< 7 , 

where M!^ is the space of matrices A, formed by 2 x 2-blocs such that 
\A\^ :=sup(a)( 6 )|A|^max([a- 6 ],iPe'^i[“-''l < 00 . 

a,b 

It is easy to see that C if = 1 and m*, entering the definition of 

is sufficiently big. So R2 also holds. 

Let us set uo{9,x) = u{0,9;0){x). Then for every I G R" and 9o G T'’^* the 
function (t,x) i->- uo(9o + tbj,x) is a solution of (17.11) with e = 0. Application of 
Theorem l3.6l and Corollary 13.71 gives us the following result: 

Theorem 7.1. For e sufficiently small there is a Borel subset 'De G T>, meas(21 \ 
Ttg) < C'e“, a > 0, such that for p G Fg there is a function ui{9,x), analytic in 
9 G T'l and H‘^ -smooth in x G , satisfying 

sup ||ui(6»,-) - uo(6»,-)|| I3e, 

|ae|<f 

and there is a mapping uj' : Vg ^ R", ||u;' — a;||ci(x)^) < fie, such that for p G Vg 
the function u{t,x) = ui{9tui'{p), x) is a solution of the beam eguation (17.11) . 
Eguation (ESI), linearised around its solution ip{t), corresponding to the solution 
u{t,x) above, has exactly N unstable and N stable directions. 

The last assertion of this theorem follows from the item (iii) of Theorem l3.6l which 
implies that the linearised equation, in the directions, corresponding to C, reduces 
to a linear equation with a coefficient matrix which can be written as B = BjrtBBao. 
The operator Bjr is close to the Hamiltonian operator with the matrix H , so it has 
N stable and N unstable directions, while the matrix Boo is skew-symmetric, so it 
has imaginary spectrum. 

Remark 7.2. This result was proved by Geng and You [9] for the case when the 
perturbation g does not depend on x and the unperturbed linear equation is stable. 

7.2. NLS equation with a smoothing nonlinearity. Consider the NLS equa¬ 
tion with the Hamiltonian 

g{u) = \ J |Vup da;-I-^ J \u{x)\'^ dxe j f{t,{—A)~°‘u{x),x)dx, 

where m > 0, a > 0 , u{x) is a complex function on the torus T'’** and / is a real- 
analytic function on R x R^ x (here we regard C as R^). The corresponding 
Hamiltonian equation is 

(7.3) ii = i{ - A-\- mu £(-A)““V 2 /(t, (-A)““m(x), x)) , 

where V 2 is the gradient with respect to the second variable, u G R^. We have to 
introduce in this equation a vector-parameter p G R”. To do this we can either 
assume that / is time-independent and add a convolution-potential term V{x, p)*u 
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(cf. (j7.1ll 'l. or assume that / is a quasiperiodic function of time, / = F{pt, u(x), x), 
where p €1) R". Cf. [2]. 

Let us discuss the second option. In this case the non-autonomous equation 
(17.31) can be written as an autonomous system on the extended phase-space O x 
T” X L2 = {(^j^(•))}! where L 2 = L2(Tr‘^*;R^) and is a ball in R”, with the 
Hamiltonian 

g{r,u,p) = hup{r,u,p)+£ j F{9,{-A)~‘^u{x),x) dx, 

hup{r,u,p) = {p,r) + ^ J \Vu\'^ dx+^ J \u{x)\'^ dx. 

Assume that m > cQ and take for A^p the operator — A -|- to with the eigenvalues 
Aa = |ap + Tn. Then the Hamiltonian g{r^ u, p) has the form, required by Theorem 
13.61 with 

/: = Z‘^*, J'=0, x = min(2a,l), ^1 = 2, /32 = 0, h=2 

(any will do here in fact) and suitable tr,/r,71 > 0 and 72 = to*. The theorem 
applies and implies that, for a typical p, equation (1731) has time-quasiperiodic 
solutions of order e. The equation, linearised about these solutions, reduces to 
constant coefficients and all its Lyapunov exponents are zero. 

If a = 0, equations (17.31) become significantly more complicated. Still the asser¬ 
tions above remain true since they follow from the KAM-theorem in [7]. Cf. [6], 
where is considered nonautonomous linear Schrodinger equation, which is equa¬ 
tion (1731) with the perturbation £(—A) “V 2/ replaced by £V{pt,x)u, and it is 
proved that this equation reduces to an autonomous equation by means of a time- 
quasiperiodic linear change of variable u. In [5] equation (1731) with a = 0 and 
/ = F{pt,{—A)~°‘u{x),x) is considered for the case when the constant-potential 
term mu is replaced by V(x)u with arbitrary sufficiently smooth potential V(x). It 
is proved that for a typical p the equation has small time-quasiperiodic solutions, 
but not that the linearised equations are reducible to constant coefficients. 

Appendix A. 


A.0.1. Transversality. 

Let V be the unit ball in R^*. For any matrix-valued function 

/ : 22 5Z(dim,C), 
let 

S(/,e) = {pG22: ||/(p)-i|| >1}, 

II II g. 

where || || is the operator norm. 

Lemma A.l. Let / : 22 —>■ C fee a -function which is j, So)-transverse, 1 < 
j < s*. 

Then, 

Leb{p G 22 : |/(p)| < e} < 

00 00 

C is a constant that only depends on s, and p. 

if undesirable, the term imu can be removed from eq. 117.3ll by means of the substitution 
u{t, x) = u'{t, 
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Proof. It is enough to prove this for 3 = (1,0,... ,0), i.e. for a scalar p. It is a 
well-known result, see for example Lemma B.l in [4], that 


Leb(S(/,e)) <C 


\fl 




This implies the claim. 

□ 


A.0.2. Extension. 

Lemma A.2. Let X C Y be subsets oJVq such that 

distfPn \ Y.X) >e, 

then there exists a C°°-function g : Eq —>■ K, being = 1 on X and = 0 outside Y 
and such that for all j > 0 

igicpvo) < 


C is an absolute constant. 

Proof. This is a classical result obtained by convoluting the characteristic function 
of X with a C^-approximation of the Dirac-delta supported in a ball of radius 
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